AO-AO«9  S65 


UNCUASSIFICO 


R AND  0 ASSOCIATES  MARINA  DEL  RCY  CALIF  F/«  19/5 

REANALYSIS  AND  APPLICATION  COMPUTER  PR09RAMS  FOR  IMPR0VIN9  ARTI— ETC(U) 
JAN  76  R L MANCUSO*  R 6 HAOFIELO  DAAB07-7A-C-0161 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BURLAIi  Of  STANDARDS  1963-^ 


Final  Report 


January  1976 


REANALYSIS  AND  APPLICATION  COMPUTER 
PROGRAMS  FOR  IMPROVING 
ARTILLERY  ACCURACY 


R.  L.  MANCUSO 


R.  G.  HADFIELD 


Prepared  for 


ATMOSPHERIC  SCIENCES  LABORATORY 
U.S.  ARMY  ELECTRONICS  COMMAND 
WHITE  SANDS  MISSILE  RANGE,  NEW  MEXICO  88002 


Contract  DAAB07-74-C-0181 
SRI  Project  3395 


D D C 

?f7apnn  fir 

FEB  3 1978 


STANFORD  RESEARCH  INSTITUTE  “ 

Menlo  Park,  California  94025  • u.S.A. 

sItatement  X 

S AppfOTW]  for  public  releoM;  I 
iDUMbutlon  Unllmitad  I 


I 


^03 


STANFORD  RESEARCH  INSTITUTE 

Menlo  Park,  California  94025  U S. A. 


Final  ^e/Srni,^ 


Jan^^  it76 


( 6 / JlEAN  ALYsis  AND  ^APPLICATION  COMPUTER 
r ^ ^ /ROGRAMS""fOR JMPROvTnG 

ARTILLERY  ACCURACY. 


,r.j,./mancuso  i 
/ .5^J^adfield  j 


Prepared  for 

ATMOSPHERIC  SCIENCES  LABORATORY 
U.S.  ARMY  ELECTRONICS  COMMAND 
WHITE  SANDS  MISSILi_fiANGE,  NEW  MEXICO  88002 


Contract  DAABO/L 74-0^1 81 
SRI  Project  3395 


^ 


' '7 


BUTTON  STAl 


Approved  ior  public  release; 
iDisMbuttos  Unlimited 


m 


r 


^ ABSTRACT 

* ^ Two  programs  have  been  developed  for  using  both  Air  Force  Global 

Weather  Central  (GWC)  prognostic  data  and  Army  rawinsonde  data  for  improving 
the  accuracy  of  meteorological  corrections  applied  in  aiming  cannon 
artillery.  The  first  program,  the  Prognostic  Data  Reanalysis  Routine 
(PDRR),  updates  the  GWC  prognoses  using  the  available  Army  rawinsonde 
observations.  The  program  options  include  data  editing,  updating  for 
both  current  and  later  times,  divergence  alteration,  and  height  balancing. 
The  second  program,  the  Artillery  Applications  Routine  (AAR),  uses  the 
updated  GWC  results  in  addition  to  the  best  available  sounding,  to  generate 
a Computer  Meteorological  Message  for  use  in  computing  ballistic  corrections 
for  artillery  firings.  The  options  that  are  to  be  used  in  either  program 
are  determined  by  control  parameter  values  that  are  initially  read  from 
data  cards  into  the  computer.  A description  of  the  numerical  techniques 
and  the  computer  program  input  and  output  are  given  in  this  report,  along 
with  a listing  of  the  programs  and  an  example  case. 
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I INTRCroUCTION 


Artillery  accuracy  is  significantly  dependent  on  accurate  and  detailed 
information  regarding  weather  conditions  (Shinn  and  Maynard,  1974).  In 
particular,  wind,  temperature,  and  density  affect  the  motion  of  a projectile 
through  the  atmosphere.  Meteorological  corrections  are  made  either  manually 
using  a Ballistic  Message,  or  by  preparing  a Computer  Met  (meteorological) 
Message  (U.S.  Army  Manual  6-15,  1970).  In  the  second  method,  a computer 
is  used  to  determine  the  ballistic  trajectories  and  the  required  meteoro- 
logical corrections.  Meteorological  information  is  normally  obtained  from 
one  dedicated  GMD-1  rawinsonde  section  that  is  assigned  to  a division  or  corps 
artillery.  Various  new  methods  are  currently  being  considered  by  the 
U.S  Army  for  improving  the  accuracy  of  meteorological  information  by  basing 
it  on  five  or  more  rawinsonde  observations  (Barnett  et  al.,  1974).  These 
methods  include: 

• Best  available — This  method  is  based  on  using  the  observation 
with  the  minimum  time-distance  deterioration  factor.  The 
variability  over  30  km  is  approximately  equivalent  to  that  for 
a 1-hour-old  message. 

• Weighted  mean--In  this  method  the  weighting  is  based  on  the  time 
and  distance  of  observation  from  the  artillery  firing.  This 
method  even  with  equal  weighting  proved  superior  to  the  conven- 
tional approach  in  tests  based  on  both  sounding  and  artillery  data 
of  the  NATO  "Summerwind"  exercise  at  Meppen,  Germany. 

• Extrapolated/ Interpolated — This  method  includes  several  approaches, 
among  them  the  use  of  a third-degree  polynomial  spline  for  extra- 
polation of  meteorological  conditions  forward  in  time  and  the  use 
of  a smoothed  plane  for  interpolation  in  area. 

• GWC/ARTY  forecast--This  method  uses  the  USAF  Global  Weather  Central 
(GWC)  numerical  prognostic  products  (U.S.  Army  Force,  1973),  and 
updates  them  with  local  U.S.  Army  artillery  rawinsonde  data.  It 

is  the  method  described  in  this  report. 


★ 


References  are  listed  at  the  end  of  the  report. 


The  GWC/ARTY  forecast  method  provides  a very  general  and  versatile 
approach,  since,  by  suitable  adjustment  of  parameters,  it  can  give  results 
equivalent  to  those  of  the  other  methods.  However,  by  using  the  GWC 
products,  more  reasonable  meteorological  information  should  be  obtained  for 
future  times  and  for  data-void  areas.  Accurate  or  just  reasonable  infor- 
mation regarding  the  weather  over  enemy  territory  would  be  critical  in 
actual  battle  situations.  Preliminary  testing  of  this  forecast  method 
was  carried  out,  using  both  rawinsonde  observations  of  the  National  Severe 
Storms  Laboratory  (NSSL)  and  of  the  U.S.  Army  Artillery  Met  Comparisons. 
The  Artillery  Met  Comparisons  were  made  at  White  Sands  Missile  Range  (WSMR) 
by  the  Atmospheric  Sciences  Laboratory  (ASL)  of  the  U.S.  Army  Electronics 
Command,  during  November  and  December  of  1974  (ASL,  1974).  This  report 
presents  results  obtained  for  one  of  these  days--19  November  1974.  Results 
for  the  test  based  on  the  NSSL  data  have  been  given  in  Quarterly  Report  3 
(Mancuso,  1975). 


Special  upper-air  observations  have  been  carried  out  by  NSSL  since  1966. 
The  upper-air  mesonetwork  for  1966  and  1967  consisted  of  11  sites 
separated  by  about  50  miles.  Soundings  were  generally  made  up  to  100  mb 
and  at  90-minute  Intervals  during  periods  of  expected  severe  weather. 
Data  for  these  years  were  obtained  from  NSSL  (Barnes  et  al.,  1971). 


II  PROGNOSTIC  DATA  REANALYSIS  ROUTINE  (PDRR)  j 

The  PDRR  program  is  for  reanalyzing  or  updating  the  QJC  prognostic 
data,  using  the  most  current  Army  artillery  upper-air  observations.  The  | 

data  treated  are  the  zonal  and  meridional  wind  components  (U  and  V),  virtual 
temperature  (T) , and  pressure  height  (H) . Figure  1 shows  the  standard 
GWC  polarstereographic  grid  for  the  northern  hemisphere.  Over  North  America, 

Europe,  and  Asia,  mesoscale  grids  are  used  with  mesh  lengths  half  the 
standard.  The  PDRR  program  uses  a 5-by-5  section  of  either  the  Standard 
or  mesoscale  grid.  The  C^C  grid  section  over  White  Sands  used  in  this 
study  is  indicated  in  both  Figures  1 and  2.  Since  this  lies  within  the 

North  American  mesoscale  grid,  it  covers  only  a 2-l/2-by-2-l/2  section  1 

of  the  standard  grid,  and  has  grid  points  about  80  nautical  miles  apart. 

When  a set  of  the  GWC  prognostic  data  becomes  available,  the  motions 

I 

or  trends  of  the  weather  patterns  are  determined  for  different  forecast 
times  and  for  the  different  GWC  prognostic  levels  (surface,  850,  700,  and 
500  mb) . The  tracking  technique  used  to  accomplish  this  was  developed 
by  Wolf  (Mancuso  and  Wolf,  1974).  This  technique  consists  of  selecting  a 
3-by-3  matrix  of  data  values  (template)  from  the  QUC  5-by-5  grid  of  data 
for  a given  time.  The  template  is  then  placed  at  some  possible  location 
in  the  data  at  the  next  later  time,  and  in  this  application  an  rms  difference 
is  computed  between  the  template  and  the  corresponding  template  location 
for  the  later  time.  The  search  for  a minimum  rms  starts  with  the  template 
in  its  initial  position,  assuming  zero  motion;  an  rms  difference  is 
computed  and  used  as  the  local  minimum.  Then  rms  values  are  computed 
for  the  template  locations  to  the  right,  above,  below,  and  the  left. 
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When  a lower  rms  is  found,  the  minimum  is  moved  to  the  location  with  the 
lower  rms,  and  the  search  begins  at  the  new  minimum.  When  this  search 
ends,  with  no  new  minimum  having  been  found,  all  locations  within  one 
element  of  the  final  minimum  are  checked.  If,  in  this  process,  a lower 
rms  value  is  encountered,  the  minimum  is  moved,  and  the  search  is  started 
at  the  new  location.  In  this  manner,  all  locations  in  a neighborhood  of 
the  minimum  are  checked,  and  the  rms  value  is  certain  to  be  a local  minimum. 
The  result  of  applying  the  above  technique  is  a motion  in  whole  grid  units. 
Thus,  to  get  a more  accurate  motion,  it  is  necessary  to  interpolate.  The 
interpolation  method  used  is  an  application  of  a two-dimensional  parabola 
in  both  directions.  The  difference  between  the  initial  location  of  the 
template  and  the  final  interpolated  location  gives  a motion  or  trend  in 
grid  units,  which  can  be  converted  to  meters  per  second  by  knowing  the 
latitude  and  longitude  of  the  initial  template  location,  the  grid  spacing, 
and  the  time  between  data  sets. 

In  the  next  step,  the  GWC  prognostic  data  are  interpolated  onto  a 
relatively  small  spherical  grid  that  covers  the  area  of  military  interest. 
This  minigrid  is  used  for  performing  all  the  remaining  computations.  An 
example  of  a minlgrid  for  the  White  Sands  area  is  shown  in  both  Figures 
2 and  3.  The  actual  size  of  the  minigrid  is  controlled  by  input  parameters; 
those  chosen  would  depend  on  the  actual  military  situation. 

The  Army  observations  that  are  available  for  analysis  are  first 
edited  for  errors  or  nonrepresentative  wind  values  by  comparing  each  datum 
with  adjacent  data,  including  the  WC  prognostic  values.  The  editing  is 
based  on  a method  developed  by  Endlich  et  al.  (1972).  It  is  accomplished 
automatically  be  computing  a wind  vector  at  the  location  of  the  measure- 
ment, using  the  technique  used  for  computing  a grid-point  wind  (Appendix  A). 
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The  editing  is  based  principally  on  surrounding  wind  data,  however,  the 
GWC  grid-point  wind  value  nearest  the  location  is  included  in  the 


computation  at  a relatively  low  weight.  The  analyzed  wind  vector  Wg 
is  then  compared  with  the  measured  wind  vector  by  taking  the  dot 
product:  e = * IVa/lV^max  (of  a and  m) . This  provides  a convenient  form 

for  comparing  vectors.  The  computation  is  similar  to  that  for  a vector 
correlation  coefficient;  when  the  two  vectors  are  identical,  e is  equal 
to  one,  but  when  they  differ,  e is  less  than  one.  If  e is  less  than  some 
preset  value  e^,  the  measurement  is  considered  inconsistent  with  other 
data  and  is  deleted.  It  is  advantageous  to  set  the  criterion  so  that 
somewhat  more  data  are  deleted  than  merely  the  bad  data.  Although  this 
tends  to  omit  good  values,  it  ensures  the  elimination  of  almost  all  the 
bad  ones. 

The  Army  observations,  which  are  grouped  within  a specified  time 
interval,  are  used  to  update  the  GWC  prognostic  fields  for  the  midpoint 
of  the  time  interval.  The  data  are  assumed  to  be  representative  of  what 
actually  occurs  over  the  release  site  at  the  initial  time  of  the  balloon 
launching.  This  approach  was  used,  since  it  is  as  reasonable  as  attempting 
to  correct  for  balloon  displacements  and  having  the  analysis  for  different 
levels  apply  to  different  times  (see  Appendix  B) . Two  methods  for  u(>dating 
the  GWC  prognostic  fields  were  developed  and  tested. 

In  the  first  method  (Update  1),  a grid-point  or  objective  analysis 
is  made,  using  both  the  U.S.  Army  observations  and  the  original  GWC  fields. 
The  new  or  updated  grid-point  values  are  computed  by  a least- squares 
fitting  of  a first-degree  polynomial  to  nearby  data.  In  the  least-square 
fitting,  the  observations  closest  to  the  grid  point  are  given  the  greatest 
weighting,  and  a greater  weighting  is  also  given  to  up-  and  down-stream 
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observations  than  to  cross-stream  ones;  this  procedure  produces  an 
anisotropic  or  elliptical-shaped  weighting  pattern.  The  initial  GWC 
grid-point  value  is  also  included.  It  is  assigned  a relatively  low  weight, 
so  that  it  will  be  significant  only  in  areas  with  no  nearby  observations. 

This  analysis  method  has  generally  given  very  satisfactory  results;  it  has 
been  described  in  detail  by  Mancuso  and  Wolf  (1974')  and  is  also  described 
in  Appendix  A. 

In  the  second  method  (Update  2),  the  GWC  minigrid  fields  of  U,  V,  H, 
and  T are  moved  an  equal  amount  until  they  give  a best  fit  to  the  observed 
Army  observations.  The  best  fit  occurs  when  the  normalized  mean-squared 
error  for  all  the  fields  combined  has  been  minimized.  The  grid  movements 
are  done  independently  for  each  level.  This  method,  when  tested  using 
the  NSSL  data,  showed  large  Inconsistencies  in  the  results  between  different 
levels.  Thus,  further  testing  and  improvement  are  required  before  applying 
this  update  method,  and  it  is  not  currently  included  in  the  program  options. 

When  the  GWC  fields  have  been  updated,  a wind  field  is  computed  that 
retains  the  updated  vorticlty  field  but  has  the  original  GWC  divergence  values, 
using  the  wind-altering  technique  of  Endlich  (1967) . The  vorticity  field  is 
computed  from  the  updated  grid-point  values  of  the  U and  V wind  components. 

The  altered  winds  are  then  used  to  compute  the  balance  equation  term  at 
each  grid  point,  and  these  terms  are  used  to  compute  a balance-height 
field.  Normal-gradient-boundary  conditions  and  an  alternating-direction- 
implicit  relaxation  method  are  used  to  perform  this  computation  of  balanced 
heights  (Mancuso,  1967).  Results  for  both  the  original  and  the  final 
updated  GWC  fields  are  discussed  in  Section  IV-C. 
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Difference  fields  between  the  original  and  the  updated  GWC  prognostic 
fields  are  computed  for  each  level.  These  difference  fields  are  moved 
with  the  GWC  pattern  trends,  using  the  upwind  advection  scheme.  The 
advected  difference  fields  are  then  used  to  adjust  or  update  the  GWC 
prognoses  for  later  times.  Thus,  the  basic  steps  of  the  PDRR  computer 
program  are  as  follows; 

• When  a new  set  of  CMC  analysis  and  prognostic  data  (U,  V,  T,  and 
H)  first  becomes  available,  it  is  used  to  determine  the  motion 
and  trends  of  the  weather  patterns.  This  is  done  by  using  a 
pattern- tracking  technique. 

• The  GWC  data  are  interpolated  onto  a minigrid  that  is  used  for 
performing  all  the  remaining  computations.  This  minigrid  is  a 
spherical  coordinate  grid  system  that  covers  the  area  of  military 
interest. 

• U.S.  Army  rawinsonde  data  (U,  V,  T,  and  H)  are  used  to  update  the 
OJC  prognostic  data.  The  Army  observations  are  grouped  within  spec- 
ified time  Intervals,  and  are  treated  as  sets  of  data  occurring  at 
the  midpoint  of  the  time  intervals.  The  observed  wind  data  are 
checked  for  errors  by  comparing  each  datum  with  adjacent  obser- 
vations and  the  nearest  GWC  value. 

• The  first  observational  set  of  data  at  time  t|^  is  used  to  update 
the  (SJC  prognostic  fields  for  that  same  time,  as  described 
previously.  Difference  fields  between  the  original  GWC  prognostic 
and  the  updated  fields  are  computed.  These  difference  fields  are 
moved  with  the  speed  and  direction  computed  initially  for  the 

GWC  patterns  and  are  used  to  adjust  the  prognostic  field  for  times 
later  than  t^. 

• Second  and  subsequent  sets  of  Army  observational  data  are  similarly 
used  to  further  update  the  latest  updated  GWC  fields. 

• When  a new  set  of  GWC  data  is  received,  the  above  procedure  is 
repeated.  However,  only  the  new  (^C  data  are  used,  and  all 
previous  values  are  purged  from  the  active  data  base. 


In  the  current  computer  program  the  time  interval  within  which 
the  Army  observations  are  grouped  is  set  at  40  minutes  or  ±20 
minutes  (IDIF=20) . 
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III  ARTILLERY  APPLICATIONS  ROUTINE  (AAR) 

The  purpose  of  the  AAR  program  is  to  provide  the  meteorological  in- 
formation in  the  form  required  for  artillery  firings.  The  program  uses 
the  PDRR  output  that  contains  the  GWC  updated  fields  on  the  minigrids. 
Values  of  U,  V,  H,  and  T at  the  four  GWC  levels  (surface,  850,  700,  and 
500  mb)  are  first  interpolated  from  the  GWC  updated  data  for  the  specific 
time  and  location  of  the  firing.  A more  reasonable  profile  for  each  of 
the  four  quantities  is  then  constructed  that  fits  the  updated  GWC  values 
at  the  four  levels,  but  that  retains  the  shape  of  the  best  available 
observed  profile.  Figure  4 shows  a typical  U component  profile  (solid 
line)  that  has  been  derived  by  shifting  the  best  available  profile 
(dashed  line)  to  fit  updated  data  that  are  indicated  by  an  asterisk  (*) ; 
the  derived  curve  goes  through  the  updated  data  points.  The  error 
quantities,  E^,  indicate  the  difference  between  the  best  available  profile 
value  and  the  updated  value.  The  solid  curve  is  obtained  by  subtracting 
the  error  value  from  the  best  available  values  at  the  WC  levels,  and 
by  subtracting  a weighted  average  error  value  for  points  in  between  the 
GWC  levels.  The  modified  profile  values  are  then  converted  into  a Computer 
Met  Message  in  terms  of  height  zone  values  for  use  in  computing  the 
artillery  corrections.  Table  1 shows  an  example  of  a Computer  Met 
Message  format  as  used  by  artillery  units. 


The  "best  available"  is  defined  as  the  observation  with  the  lowest  At, 
where  At  = (time  in  minutes  between  observation  and  firing)  plus  (distance 
in  kilometers  from  gun  to  observation,  divided  by  2) . 
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IV  COMPUTER  PROGRAMS 

A flow  diagram  for  the  H)RR  and  AAR  computer  programs  is  shown  in 
Figure  5.  The  programs  were  originally  written  for  a CDC  6400  computer, 
and  FORTRAN  listings  compatible  with  that  computer  are  given  in  Appendices 
C and  D.  The  programs  were  modified  so  that  they  could  be  run  on  the 
WSMR  Univac  1108  computer.  These  versions  of  the  programs  have  been  given 
to  the  Atmospheric  Sciences  Laboratory  of  the  U.S.  Army  Electronics 
Command  at  WSMR.  The  programs  were  designed  to  be  as  practical  as  possible 
in  regard  to  the  data  they  require  and  the  ease  with  which  they  can  be 
modified  to  incorporate  improved  products  and  the  results  of  continuing 
research.  Thus,  the  programs  were  organized  in  the  form  of  several 
subprograms,  and  changes  can  be  made  with  a minimum  of  disruption  of  the 
rest  of  the  program  . If  future  research  provides  better  methods  for 
updating  OlC  fields,  they  can  be  incorporated  relatively  easily. 

A.  Data  Input  for  the  FPRR  Program 

The  input  data  card  format  for  the  PDRR  program  is  shown  in  Table  2 
with  example  values.  The  input  data  cards  provide  a convenient  way  for 
reading  in  new  values  for  the  control  parameters  of  the  program.  The 
input  parameters  on  Card  A (ICHK,  lANA,  IDIV,  IBAL,  IFOR,  IPRINT,  and 
ITAPE)  are  control  parameters  that  specify  whether  or  not  a specific  task 
is  to  be  performed.  That  is,  if  the  parameter  value  is  set  greater  than 
0,  the  task  is  performed;  if  it  is  set  equal  to  or  less  than  0,  the  task 
is  not  performed.  Specifically,  if  ICHK  is  set  greater  than  0,  then 
subroutine  CHECK  is  called,  and  the  observed  wind  values  are  edited  by 
comparing  them  with  analyzed  values  for  the  same  locations  (if  an  incon- 
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sistency  is  found  between  an  observed  and  an  analyzed  wind,  the  observed 
wind  is  deleted) . IF  lANA  is  set  greater  than  0,  the  Army  observed  data 
are  analyzed  with  the  GWC  data,  using  Update  1.  If  IDIV  is  set  greater 
than  0,  the  divergence  values  of  the  wind  fields  are  made  equal  to  the 
aJC  values.  If  IBAL  is  set  greater  than  0,  height  fields  are  computed 
that  are  balanced  with  the  wind  fields.  If  IFOR  is  set  greater  than  0, 
forecast  fields  are  generated.  The  symbols  IPRINT  and  ITAPE  are  control 
parameters  for  the  output.  If  IPRINT  is  set  greater  than  0,  the  output 
is  printed  out  on  the  line  printer.  If  ITAPE  is  set  greater  than  0,  the 
output  is  printed  on  tape. 

Card  B contains  information  regarding  the  size  and  location  of  the 
minigrid  array.  N9  indicates  the  number  of  columns,  and  M9  indicates  the 
number  of  rows.  XB  and  YB  are  the  longitude  and  latitude  locations  of 
the  upper-left-hand  corner  of  the  minigrid,  and  the  symbol  DD  indicates 
the  distance  between  the  grid  points  in  degrees.  A different  minigrid 
is  used  in  the  computer  example  case  then  that  illustrated  in  Figures  2 
and  3.  XC  and  YC  establish  the  upper- left-hand  starting  location  of  the 
5-by-5  section  of  the  GWC  grid.  The  GWC  data  is  read  in  from  a binary 
tape  on  unit  1 in  the  subroutine  GWCIN  (Appendix  C) ; using  two  separate 
read  statements. 


Input  data  are  given  on  Card  C for  the  symbol  IT,  which  represents 
the  total  number  of  Army  data  sets  to  be  analyzed.  Cards  D contain  the 
input  information  for  array  LDATE,  that  is,  the  data  and  midtime  for 
each  of  the  Army  data  sets  to  be  used  for  updating.  Although  values 
for  only  one  D card  are  used  in  the  example  case  (IT=1) , IT  may  be  increased 
to  any  desired  value;  each  date/ time  would  then  be  specified  on  a separate 
D card  in  the  same  format  as  shown  in  Table  2.  The  Army  upper-air  data 
of  WSMR  is  read  in  from  cards  by  the  RAWIN  subroutines  of  both  the  PDRR 
and  AAR  programs  (Appendices  C and  D),  using  a 10F8.1  format. 
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Card  E is  used  to  designate  the  number  of  Army  rawinsonde  stations 
(JT)  that  are  being  used.  Cards  F contain  input  values  for  the  locations 
and  elevations  of  the  total  number  of  stations  that  are  specified  on  the 
E card.  Each  F card  contains  data  for  one  of  the  JT  stations  for  the 
variables  XS(J),  YS(J),  and  ES(J),  each  integer  of  J representing  one 
of  the  JT  stations.  An  example  of  input  data  for  the  PDRR  program  is 
given  in  Appendix  E. 


Data  Input  for  the  AAR  Program 


A description  of  the  input  data  card  for  the  AAR  program  is  given  in 
Table  3 along  with  example  values.  Card  A contains  a list  of  the  various 
control  parameter  symbols  that  are  used  for  printing  out  labels  in  the 
program.  Corresponding  numerical  values  for  each  of  the  parameters  are 
read  in  on  Card  B.  The  parameter  IQ  is  an  indicator  that  controls  what 
final  data  are  to  be  printed  out.  If  IQ  is  1 or  greater,  profile 
curves  are  plotted  out,  and  a Computer  Met  Message  is  printed  out  in  tabulated 
form.  If  IQ  is  less  than  1,  then  only  the  Computer  Met  Message  is  printed 
out.  The  variable  XFIRE  is  the  longitudinal  location  of  the  center  of 
the  firing  path.  The  variables  XUL,  XD,  and  NX  refer  to  the  minigrid.  XUL 
is  the  longitude  of  the  upper- left-hand  comer  of  the  minigrid;  XD  is  the 
distance  in  degrees  of  longitude  between  the  grid  points;  and  NX  is  the 
F number  of  grid  points  in  the  X direction  (longitude).  The  variable  YFIRE 

' is  the  latitudinal  location  of  the  center  of  the  firing  path.  The  variables 

I YUL,  YD,  and  NY  also  refer  to  the  minigrid.  YUL  is  the  latitude  of  the 

* upper-left-hand  comer  of  the  minigrid;  YD  is  the  distance  in  degrees 

latitude  between  grid  points ; and  NY  is  the  number  of  minigrid  points  in 
; the  Y direction  (latitude) . 
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Card  C contains  a listing  of  the  control  parameter  symbols  for  the 
values  given  on  Card  D.  The  parameter  NT  of  Card  D is  the  number  of  GWC 
forecast  times,  including  the  initial  analysis  (NT=19).  NL  is  the  number 
of  levels  that  are  available  from  the  updated  Q^C  data  (NL=4) . TFCST  is 
the  time  (GMT)  of  the  initial  analysis  of  the  QJC  data,  and  TFIRE  is  the 
time  (GMT)  of  the  artillery  firing.  IDATE  specifies  the  year,  month,  and 
day  of  the  firing.  JSTAT  is  the  number  of  Army  observation  stations  that 
are  used  for  determining  the  best  available  sounding. 

Card  E contains  a listing  of  the  station  longitude,  latitude,  and 
elevation  symbols.  Cards  F contain  values  for  these  variables,  that  is, 

XS(J)  in  degrees  longitude,  YS(J)  in  degrees  latitude,  and  ES(J)  in  meters 
above  sea  level;  each  integer  of  J represent  one  of  the  JT  stations.  An 
example  input  for  the  AAR  program  is  given  in  Appendix  F. 

C . Data  Output 

The  data  output  from  both  the  PDRR  and  the  AAR  programs  can  be  controlled 
with  the  IPRINT  and  ITAPE  parameters  of  Table  2 and  the  IQ  parameter  of 
Table  3.  With  the  PDRR  program,  detailed  line  printer  outputs,  as  shown 
in  Appendix  G,  can  be  obtained  (IPRINT=1) , or  tape  output  only  can  be 
obtained  (IPRINT=0,  ITAPE=1).  In  the  print  output  of  the  PDRR  program. 

Levels  1,  2,  3,  and  4 refer  to  the  surface,  850,  700,  and  500  mb.  With 

the  AAR  program,  one  can  choose  to  obtain  the  curves  showing  the  best 
available  and  updated  profiles,  in  addition  to  a Computer  Met  Message, 

as  shown  in  Appendix  H (IQ=1).  If  desired,  however,  only  the  Computer 

Met  Message  need  be  printed  out  (IQ=0).  The  storage  requirement  of  the 
PDRR  program  is  23,000  words  of  memory,  and  the  central  processing  time 
for  the  example  used  in  Appendices  E and  G is  26  seconds  (CDC-6400).  The 
storage  requirement  of  the  AAR  program  is  16,000  words  of  memory, 
and  the  central  processing  time  for  the  example  of  Appendices  F and  H 
is  12  seconds. 
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V SUMMARY  AND  RECOMMENDATIONS 

A method  for  providing  meteorological  information  for  artillery 
has  been  developed  in  this  study  and  implemented  as  the  PDRR  and  AAR 
computer  programs.  The  first  of  these  programs  uses  artillery  upper-air 
observations  to  update  GWC  prognostic  fields.  The  second  generates 
atmospheric  zone  values  of  wind,  temperature  (virtual),  and  pressure 
for  the  Computer  Met  Message,  on  the  basis  of  updated  GWC  fields  and 
the  best  available  rawinsonde  profile.  These  programs  provide  a versa- 
tile and  general  method  for  supplying  meteorological  information  to  the 
artillery,  and  would  be  particularly  applicable  where  extrapolation  is 
required  in  either  time  or  space.  The  method  uses  the  current  numerical 
products  of  GWC,  which  are  based  on  a baroclinlc  prediction  model.  How- 
ever, forecast  should  eventually  be  based  on  a primitive  equation  (PE) 
model,  which  would  use  smaller  grid  spacing  and  more  levels.  Also,  a 
mesoscale  forecasting  model  has  been  developed  by  the  U.S.  Air  Force 
(Kaplan,  1972),  Results  from  either  of  these  would  be  more  appropriate 
to  use  as  the  basic  fields  to  be  updated  with  the  Army  observations  and 
should  lead  to  substantial  improvements.  This  would  be  particularly 
important  for  providing  more  accurate  surface  forecasts  over  terrain 
such  as  White  Sands.  If  either  of  the  above-mentioned  types  of  models 
are  made  operational  by  GWC,  the  PDRR  and  AAR  programs  should  be  modified 
to  use  the  new  products.  Other  recommendations  for  refining  the  programs 
are : 


25 


* Optimize  parameters,  and  compare  different  options  of  the 
Prognostic  Data  Reanalysis  Routine  (IDRR)  and  the  Artillery 
Application  Routine  (AAR) . 

• Develop  and  test  improved  analysis  techniques,  including 
the  incorporation  of  terrain  influence. 

The  testing  and  development  could  be  carried  out  by  using  the  entire 
set  of  the  ASL  Artillery  Meteorological  Comparisons  data  from  November 
and  December  1974.  The  actual  artillery  firing  results  that  were 
obtained  concurrently  with  the  upper-air  observations  could  be  used 
to  verify  the  meteorological  analyses.  These  refinements  of  the  programs 
could  lead  toward  a fully  automated  meteorological  service  that  would 
be  both  accurate  and  versatile,  as  required  for  military  operations 
(ASL,  1973) . Significant  progress  should  be  possible  because  of  the 
availability  ot  the  Artillery  Meteorological  Comparisons  data  and  the 
current  ASL  software  such  as  described  in  this  report. 
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Appendix  A 


OBJECTIVE  ANALYSIS 


The  fitting  of  a polynomial  surface  to  observations  adjacent  to  a 
grid  point  is  a natural  approach  for  interpolating  a grid-point  value. 
Schemes  based  on  this  approach  have  been  used  by  a number  of  investiga- 
tors, one  of  the  earliest  being  Gilchrist  and  Cressman  (1954).  In  the 
objective  analysis  described  in  this  appendix,  a grid-point  value  for  a 
quantity  is  obtained  by  a least-squares  fit  of  a first-degree  polynomial 
surface  to  adjacent  observational  values.  In  the  least-squares  fitting, 
the  observations  are  weighted  Inversely  with  distance.  Upstream  and 
downstream  observations  are  given  greater  weight  than  cross-stream 
observations.  This  anisotropic  weighting  produces  an  elliptical  field  of 
weights,  with  the  long  axis  in  the  direction  of  flow;  it  is  particularly 
desirable  when  analyziug  such  phenomena  as  jet  streams.  Sasaki  (1971) 
showed  by  the  variational  method  that  such  weighting  is  theoretically 
consistent  with  the  advection  terms  in  the  equations  of  motion.  The 
basic  objective  analysis  was  originally  developed  for  analyzing  surface 
and  upper-air  data  over  the  United  States  (Endllch  and  Mancuso,  1968). 

It  has  since  been  modified  to  allow  for  more  general  applications 
(Endllch  and  Mancuso,  1973). 

The  actual  objective  analysis  procedure  follows!  A grid-point  value 
for  a quantity  is  determined  by  fitting  a first-degree  polynomial  surface 
by  least  squares  to  K nearby  observations.  An  optimum  fit  is  obtained  by 
minimizing 
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Q = r w(q  - , 

1 

where  w is  a weighting  factor,  q is  an  observed  value  at  some  location 

(x,y),  and  q is  the  polynomial  estimate  for  the  same  location  (a  + bx  + cy). 

If  an  initial  guess  value  for  a grid  point  is  provided,  it  is  treated  as  the 

I first  observational  value  (K  * 1)  and  is  given  a fixed  weight  w^ . A rela- 

I tively  low  value  is  normally  assigned  to  w. , so  the  initial  estimate  will 

i ^ 

[ have  a significant  influence  only  in  areas  where  there  are  no  nearby  data. 

[The  remaining  K - 1 observations  are  normally  those  nearest  the  grid  point. 

However,  it  is  required,  if  possible,  that  at  least  one  observation  be 
I selected  from  each  of  the  four  angular  quadrants  (0-90,  90-180,  180-270,  and 

270-360  degrees).  This  tends  to  keep  the  analysis  continuous  in  regions  of 
poorly  distributed  data  and  large  gaps.  The  value  used  for  K has  varied 
between  3 and  10;  the  higher  value  is  used  in  areas  of  unevenly  distributed 
data. 

The  weighting  for  an  observation  is  given  by 

2 2 2 2 
w = C / [C  + R + (aS)  ] 

S = ( Ik  • IR  X W)/  I W 1 . 

The  R in  this  equation  represents  the  magnitude  of  the  horizontal  distance 
vector  (|R)  between  the  grid  point  and  the  observation.  Thus,  the  closer 
an  observation  is  to  the  grid  point,  the  greater  the  weighting.  The  degree 
of  this  effect  is  dependent  on  the  value  assigned  to  the  constant  C.  The 
S in  the  above  equation  is  Included  to  give  a greater  weighting  to  upstream 
and  downstream  observations  than  to  cross-stream  observations,  at  equal 
distances.  This  produces  the  anisotropic  weighting  pattern  mentioned 
previously.  The  control  parameter  a is  Included  so  that  the  relative 
Influence  of  the  anisotropic  weighting  can  be  varied. 
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The  result  obtained  by  fitting  a second-degree  polynomial  surface  to 
distance-weighted  U and  V component  data  at  300  mb  over  the  United  States 
is  shown  in  Figure  A-l(a),  This  figure  shows  the  isotachs  that  were  drawn 
to  the  computed  wind  speed  values  at  the  grid  points  spaced  every  2.3 
degrees  of  latitude  and  longitude  (the  results  for  wind  direction  are  not 
shown).  The  figure  also  shows  the.  basic  wind  data  in  vector  form.  The 
isotach  field  that  was  obtained  from  a direct  hand  analysis  of  wind  speeds 
is  shown  in  Figure  A-l(b).  There  are  some  noticeable  differences  between 
these  two  fields  [Figure  A-l(a)  and  (b)],  and  the  rms  difference  is  3.5  m 
sec  These  differences  are  due  in  large  part  to  the  fact  that  the  poly- 
nomial analysis  was  based  on  u v wind  component  data,  while  the  hand 
analysis  was  made  directly  from  wind  speed  data. 

Figure  A-l(c)  shows  the  isotach  field  obtained  by  fitting  a first- 
degree  polynomial  surface  to  distance-weighted  U and  V wind  component  data, 
both  including  and  excluding  the  anisotropic  weighting  effect.  A closer 
fitting  of  the  nonlinear  features  is  obtained  by  using  the  second-degree 
polynomial  [Figure  A-l(a)].  However,  the  use  of  higher-order  polynomials 
requires  more  computation  time  and  complex  programning.  Although  the 
overall  rms  difference  between  the  analyses  depicted  in  Figure  A-l(c)  is 
not  large  (1.3  m sec  ^),  the  inclusion  of  the  anisotropic  weighting  does 
produce  a desirable  streakiness  in  Jet  stream  areas  and  a better  fit  to  the 
nonlinear  features. 

Figure  A-l(d)  shows  the  isotach  field  obtained  when  a first-degree 

polynomial  surface  was  fitted  directly  to  wind  speeds,  as  opposed  to 

fitting  the  U and  V components;  in  both  cases,  optimum  weighting  values 
2 

were  used  (0=2  and  a = 2.5).  A very  good  agreement  with  the  hand  analysis 
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shown  in  Figure  A-l(b)  was  obtained  in  this  way  (rms  difference  = 3.0  m sec  ). 
However,  the  results  obtained  by  fitting  the  U and  V wind  components 
[Figure  A-l(c)]  are  in  better  agreement  with  the  second-degree  polynomial 
analysis  of  Figure  A-l(a)  (rms  difference  3.5  msec"^) , than  with  the  hand 
analysis  of  Figure  A-l(b)  (rms  difference  = 4.0  msec"^). 

The  weighting  parameters  C and  a are  determined  experimentally  to 

optimize  the  results.  A better  fit  to  nonlinear  features  occurs  as  C is 

decreased,  but  little  improvement  occurs  when  C is  reduced  excessively. 

The  distance  weighting  is  Important;  if  it  is  not  incorporated  or  if  C 

is  set  too  high,  the  analysis  is  overly  smoothed.  However,  some  smoothing 

is  desirable  because  of  data  errors  and  to  produce  a more  continuous  field. 

Thus,  the  value  of  C is  usually  not  set  too  low,  but  is  generally  set  less 

than  the  value  of  the  grid  spacing.  Figure  A-2  shows  the  rms  difference  of  a 

first-degree  polynomial  analysis  relative  to  both  a hand  and  a second-degree 

2 

pol)momlal  analysis  as  a function  of  or  (C  =2).  A definite  minimum  in  the 
differences  occurs  when  the  anisotropic  weighting  parameter  d,  is  set  at  a 
value  between  2 and  3.  Although  the  anisotropic  weighting  effect  is  not  as 
strong  as  the  distance  weighting  effect,  it  does  produce  important  differences 
in  critical  regions,  as  shown  in  Figure  A-l(c),  and  is  easily  incorporated 
into  the  analysis  scheme. 
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FIGURE  A-2  RMS  WIND  SPEED  DIFFERENCES  BETWEEN  FIRST-  AND 
SECOND-DEGREE  POLYNOMIAL  ANALYSES  (SOLID  LINE) 
AND  BETWEEN  FIRST-DEGREE  POLYNOMIAL  AND  HAND 
ANALYSES  (DASHED  LINE) 
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Appendix  B 

SIGNIFICANCE  OF  BALLOON  DISPLACEMENTS 

The  rawlnsonde  balloon  displacements  from  the  release  sites  are 
normally  not  considered  In  conventional  synoptic-scale  analyses.  The 
displacement  of  the  balloon  at  a height  H will  be  given  by  D — V t , where 
V Is  the  average  speed  between  the  surface  and  H,  and  t Is  the  time  It 
took  the  balloon  to  rise  to  H.  In  analyses  over  the  United  States,  D 
would  be  at  most  about  60  miles  at  high  altitudes  (10  km)  and  high  winds 
(100  msec  ^) , and  would  normally  be  less  than  30  miles.  However,  In 
mesoscale  analyses.  It  becomes  a more  Important  question,  since  the 
horizontal  resolution  becomes  less  than  60  miles.  The  approximate  spacing 
between  stations  for  various  mesoscale  networks  Is  given  below. 
Thunderstorm  Project  1947  (Fujlta,  1958)  5-10  miles 

NSSL  1966/1967  (Barnes  et  al,  1971)  50  miles 

WSMR  Oct/Nov  1974  (ASL,  1974)  10-20  miles 

A three-dimensional  Illustration  of  balloon  displacements  Is  shown  In 
Figure  B-1  for  the  Thunderstorm  Project.  The  displacements  are  quite 
significant,  and  for  an  accurate  analysis  they  would  need  to  be  taken 
Into  account. 

Two  approaches  are  possible  for  handling  the  balloon  displacement 
problem: 

• Correct  the  location  according  to  a calculated  displacement 
with  Increasing  height,  assuming  some  constant  balloon  rate 
of  rise  such  as  5.5  msec~^.  The  analyses  at  the  higher  levels 
would  have  to  be  applied  to  slightly  later  times. 


33 


I 


• Assume  that  a hindcast  can  be  made  by  simply  advecting  the  winds 
backward  in  time,  using  the  average  wind  value  from  the  surface 
up  to  the  level  in  question.  Thus,  no  correction  of  locations 
would  be  required;  that  is,  the  measured  wind  value  would  be 
allowed  to  apply  to  the  location  and  time  of  the  release. 

The  second  approach  is  the  easier  to  adopt,  since  it  requires  no 
corrections.  Its  merit  depends  on  the  validity  of  the  hindcast.  Vfhen 
corrections  are  most  needed,  as  in  strong  wind  conditions,  it  would 
probably  have  less  merit.  This  approach  is  currently  being  used  in  the 
PDRR  program;  however,  further  testing  should  be  carried  out  to  compare 
the  two  approaches,  using  the  WSMR  November  and  December  1974  comparison 


Appendix  C 

PDRR  PROGRAM 
LISTING 


non  *tnn  non 


?^ecediy<$ 


PROGRAM  PORR ( INPUT#  OUTPUT • T APE  1 • TAPE? . TAPE  T t TAPF4 .PUNCH) 

this  PROGRAM  IS  THE  PROGNOSTIC  DATA  RFANALYSIS  ROUTINE  (PDRR).  IT  IS 
USED  FOR  REANALYZING  OR  UPDATING  THE  GLOPAL  WEATHER  CENTER  I GWC ) 
PROGNOSTIC  DATA  USING  THE  MOST  RECENTLY  AVAILABLE  ARMY  UPPER  AlP 
OBSERVATIONS. 

DIMENSION  AND  COMMON  STATMENTS 
INTEGER  0ATA(12000) 

D I MENS  ION  U(  1 00 ) . V(  I 00  ) .H(  100)  . T<  1 00)  , VOR ( 100).OIV(100),RAL(100) 
DIMENSION  US (25)  . VS (2S  ) .HS( ?5) , TS<  25) , XS( 25) , YS(  ?S) ,FS<  25) 
DIMENSION  X(10).Y(10).OSX(100).OSY(100) .LDATF(5 ) , IM0(  I 2) 
DIMENSION  CM(I0).TM(  10  ).TRU(  100).TpV(  100) 

DIMENSION  OATI(*O0).GATl(I00)»GAT?{t00).VAT(25).WAT(25)»GAS(I00.4) 
0 I MENS  ION  UR( I 00) . VR(  1 00 ) . HR ( I 00 ) . TR ( 1 0 0 ) 

COMMON/CIS/  IS(tOOO) 

COMMON/CCHK/  CR T. S I M .G I M , I WND 
COMMON/CGO/  M9 .N9 , I 9, MB, N8. YB. XB, DD 
COMMON/CPTS/  KS,Wl ,C 2 . RM AX , KSS5 . lOS.KSW.ALPH 

EQUIVALENCE  (U.GAS  t 1 , 1 ) ) . ( V, GASI  1,2)),<H,GAS(1.3)),<T,GAS(1.A)), 

2 (UR, DAT I ( 1 ) ) , ( VR,DAT1  ( 101  ) ) , (HR, DAT  1(201)). 

3 ( TR.OATl  ( 301 ) ) , ( VOR.OSX)  , { DI V.OSY) 

SET  BASIC  CONSTANTS  (INITIALLY  SET  CORE  TO  ZERO) 

DATA  lOIF.IGMT  /20,-A20/ 

DATA  10,11,12/0,1,?/ 

DATA  CRT,SIM,GIMS, IWNO/0, 3,0.1 ,0.05, 1 / 

DATA  NVARP,NSIZE  /2,5/ 

DATA  Wl, C2, RMAX, KS. KSS5,KSW, IDS, ALPH/ 0.10. 0.02, 10. 0,10. 5,0, 1,2.0/ 
DATA  I 20 .LI  .L4  ,LI 00 ,K4 .N5 , N25  / 1 9 , I , 4 , 1 00 . 4 , 5 , 25/ 

DATA  XP, YP.ROT.DTS, ACR/24, 0 , 26. 0, BO. 0,1. 0,0. 01 74533/ 

DATA  I MO/0,31  ,59, 90, 120, 151. 181. 212, 243, 2 73,  104,  3 34/ 

1 FORMAT!  1 OX, 711 0) 

2 FORMAT(10X,2I10,5F10.2 ) 

3 FORMAT( lOX.I 10) 

4 FORMAT ( 14X, 312, 6X,  212 ) 

5 FORMAT!  1 OX,  HO) 

6 FORMAT! lOX, 3F10. 2) 

7 F0RMAT(2X.6I 10 ,2E15.3) 

READ  IN  CONTROL  PARAMETERS 

RFAD  I,  ICHK, I ANA, lOt V.IBAL, IFOR, IPRINT, ITAPE 
READ  2,N9, M9, X9, YB, DO.XC.YC 
READ  3, IT 

READ  4, ( ( LOATE(L.  I),L=1,5),I  = 1,IT) 

READ  5,JT 

READ  6.  ( XS( J)  , YS(  J ) ,ES( J ) , J=1 , JT) 

PRINT  OUT  CONTROL  PARAMETERS 

PRINT  9000 
PRINT  9100 
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PRINT  9000 
PRINT  9050 

PR  I NT  900 1 t ICHK.  I ANA.  I OIV  t IBAL . IFOR.  IPR INT.  I TAPF 
PRINT  900? .N9.M9.XB.VB.D0.XC.YC 
PRINT  9001.  I T 
PRINT  9054 

PRINT  9004.  (ILDATEIL.  I ).L  = 1.5)«I=1  .IT) 

PRINT  9005, JT 
PRINT  9056 

PRINT  90  06,<J,XS<J),Y5(J),FS<J)  .J=1.JT) 

PERFORM  BASIC  COMPUTATIONS 

IF  ( I T.GT.25.0R. JT,GT,?5)  GO  TO  999 
IF  (M9.GT. 10.OR.N9.GT.  10)  GO  TO  999 
GI M=GI MS 
IH?0=(  I 20-1 ) *60 
XI 00=DTS*60.0 
I 9=m9*N9 
MB=M9- 1 
N6=N9-1 
N49=I 9 
N98=N49+N49 
N4  7 = N9a+N49 
K,5-K4-  1 
I0=L4*K3*N49 
no  20  N=1,N9 
?0  X( N)=XB*OD*(N-l  ) 

XAVF^ACR* (-X(N5)-RnT) 

SAVF=SIN( XAVF) 
cAve=cos(xAve) 
on  22  M=1,M9 
Y<  M)=YO-00*<M-l  ) 

ANG=ACR*Y(M) 

CM (M ) = CaS<  ANG ) 

?2  TM (M )=TAN( ANG) 

XN5=-X ( N5 ) 

CGD= I .O* Y(N5) /90. 0 
1=0 

no  24  M=1 ,M9 

TL  A-taN(  (90.0-Y(M)  ) *ACR*0 . 5 ) *3  1 . 2 042 
DO  24  N=1,N9 
I = I*l 

XLO=ACR*(-X(N)-RnT ) 

DSX( I ) = ( XP-TLA*SIN{ XLO  )-XC ) *2 « 0*1 . 0 
24  OSY( 1 )=( YP+TLA*CnS( XUO ) - YC ) * 2. 0* 1 . 0 

LOOP  THROUGH  ALL  GWC  FORECASTING  TIMES 

no  50  IG=1,I20 

READ  IN  GWC  DATA 
C 

CALL  GWCINI  IG.  I GWC , JO  A TF , J T I ME , DA T 1 , IMD) 

DO  45  L=LI,L4 
N1={L-1>*LI00 
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IP=(  ( IG-t  )*L4*L-J 
JP  = 0 

DO  4?  K=I .K4 
IF  (K.F0.K4)  IP=1P-N4P 
DO  36  N=l ,N25 
Nl=Nt4 I 
D4TN=DAT 1 ( N1  ) 

MATIN)  ■=DATN 
IF  (K.GT.2)  GO  TO  3=i 
N3  = N1 

IF  (K.EQ.2)  N3=N1-N25 
OATN=OAT I (N3 ) 

VATN  = OATl  <N3't-N25) 

watin)=oatn 

IF  (K.FQ,2)  MAT{N)=VATN 
CALL  UVCONVIDATN, VATN, -XI N5) ) 

IF  IK.F0.2)  DATN=VATN 

35  JP=JP'*-l 

GATI I JP ) = GASI JP.L  ) 

GASI JP»L  >=DATN 
6AT2I JP) =nATN 

36  VATtN)=nATN 

OO  40  1=1.19 

DL  X=D5XI  I ) 

OLY=DSVI I ) 

NX=INTIOLX> 

MY= INT I OLY ) 

NM=I MY-1 ) •NS+NX 

DLX=DLX-NX 

0LY=0LY-MY 

C 

C INTERPOLATE  TO  MINI  GRID 
C 

CALL  I NT PT I NM.N5. VAT. DLX.OLY.BAT I 

IP=TP4.1 

IBAT=1NTIBAT*1 0.0) 

IF  (K.E0.3)  IflAT= IBAT*  10000 
IFIK.NE.4)  GO  TO  40 
IFIDATAI IP) .LT.O)  IBAT=-ibaT 
IDAT  = OATAI  IP)-HBAT 
40  DATA! IP)=IBAT 
42  CONTINUE 
C 

C COMPUTE  TRENDS 
C 

IF  I IFOR.LT.1.0R.IG.LT.2)  GO  TO  45 
L IG=I  IG- 1 )*L44L 

nvars-nvarp 

IF  IL.EQ.I)  nvars=nvarp-i 

CALL  TRENOINVAPS.NSITE.TPUIL IG) .TRV(LIG).GAT1.GAT2. IG) 

TRUILIG) =TRUILIG) MCGO 

TRV<LIG)=-TRV<L1G)*CG0 

CALL  UVCONV ITRUILIG) .TRVIL IG ). XN5 ) 

TRUILIG)=TRUtLIG)/OD 
TRV(LIG)=TRVIL IG)/DO 

IF  lABSITRUILlG) ) .GT.5.0)  T RU IL I G ) = S I GN I 5 . 0 . TRU ( L I G ) ) 


[ 
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IF  (ABStTRVILIG)  ) •r.T.'S.O)  TR  V ( L I G ) = S I GN<  S . 0 , TR  V < L I G > ) 

4fi  CONTINUE 
SO  CONTINU*^ 

ENO  OF  GWC  LOOP 

PRINT  2.JP.IP 
IF  (lANA.LT.l)  GO  TO  101 

LOOP  THROUGH  ARMY  RAWINSONOF  OBSERVATION  TIMES 

00  too  IR=I , IT 

1 = IR 

L02=LDATE<2, I » 

I0ATE=LDATE(5.I  ) + 60*(LDATF(4»I  ) +2 A ♦( LD AT E( 3 . I >*IM0(L02) 
2 ♦36S*LOATEC 1, I )-366) ) 

PRINT  1 . IGWC. IDATE 
IF  ( IDATE.lt. IGWC  ) GO  TO  Ooq 
IF  ( lOATF .GT. IGWC+ IH20 ) GO  TO  RPR 
ICOMP=IGWC*X100 
DO  60  IG--2,I20 

IF  ( IDATE-ICOMP)  Se.60,60 
SB  FG2=( ICOMP-IDATE )*0. 1 /XlOO 
FG1=0. 1-FG2 
GO  TO  61 

60  ICOMP=  ICOMP'fXI  00 
GO  TO  9R9 

61  I G= I G- 1 

IL  = ( IG-l  )*!0 

PRINT  1 .ICOMP. IDATE, IG 

PRINT  6,FG1,FG2 

READ  IN  RAWINSONDE  DATA 

CALL  RAWINl JT.IR, ES.UR.VR.HR.tr, IMO. IDATE, lOIF.  IGMT) 
LOOP  THROUGH  MANDATORY  LEVELS 


DO  95  L=Ll,L4 
LJ=(L-1 )*JT 
DO  68  J=1.JT 
US IJ>=URIJ*LJ) 

VSI  J)=VR(  J+LJ) 

HS( J»=HR( J+LJ) 

TS( J) =TR( J+CJI 
N1=IL+IL-1) *N49*K3+l 
N2=Nl+N49-l 
1 = 1 

DO  70  N=NI,N2 

U( I )=DATAC  N) ♦FGl ♦DAT A{ N+I D) *FG2 
V( I )=DATA1N+N49)*FG1+0ATA( N+N49+1D)*FG2 
IHI=INT(DATA1N+N9fl  )*0.0001) 
IH2=INTCr>ATA(N  + N98+I  D)  40.0001  ) 

IT  1 = DATAI N+N9B  )-IHl4l0000 
I T2*OATA<  N+N9B  + ID>-IH241  0000 
TT1=IABS( ITI) 


U2 


ir2=lABS( I T2) 

H(I)=  IH1*FG14- IH2AF02 

T(I)=  ITi*FGl«-IT2*'FG2 

70  I=I+1 
PHINT  9000 
PRINT  9059. L 

IF  ( ICHK.LE.O.OR.L.&Q.l ) GO  TO  170 

CALL  MESMI  I 0 . J T , YS. X S . VS . US . HS . T S . Y . X , V , 0 , H , T ) 

CALL  CHECK! JT.YS.XS.VS.US. Y. X.V.U) 
170  IF  I IPR INT.LT. 1 ) GO  TO  71 
PRINT  9025 
PRINT  9057 

PRINT  9007. IJ.HSI J).TS<J).USIJ».VSCJ).J=1,JTI 
PRINT  9000 
PRINT  9060 
PRINT  9062 

PRINT  901 0. IH( I ) . I=t . 19) 

PRINT  9064 

PRINT  9010. ( T( I ). 1=1 . 19) 

PRINT  9066 

PRINT  9010. (U( I ) . 1=1 . 19) 

PRINT  9068 

PRINT  901 0. ( V( I ) « 1=1 . 19) 

71  CONTINUE 

IF  (IANA-1)  85.73.72 

72  CONTINUE 
GO  TO  85 

73  OU  173  1=1.19 

VOR( I ) = 0.0 
DIV< 1 )=0.0 

173  aAL(I)=0.0 

CALL  K IDI  10. I 1 . 10.VOR.OIV.8AL.U.V.CM. TM) 

CALL  MESH! lO.JT.YS.XS.VS.US.HS.TS.Y.X.V.U.H.T) 

CALL  MESH!I2.JT.YS.XS.VS.US.HS.TS.Y.X«V.U.H.T) 

IF  (lOIV.LT.l)  GO  TO  75 
C 

C ADJUSTMENT  OF  WIND  FIELD  TO  THE  GWC  DIVERGENCE  FIELD 
C 

CALL  KID!  I 1 . I 1 . lO.VUR.BAL.OIV.U.V.CM.TM) 

CALL  ALTERSI25.0.0.0.5.TM.  CM.  U.V.VOR.DIV) 

75  IF  IIBAL.LT.l)  GO  TO  85 
C 

C COMPUTATION  OF  BALANCED  HEIGHT  FIELD 
C 

CALL  KIDIll.Il.il . VOR.OIV .BAL.U.V.CM. TM) 

SUH=0.0 
DO  80  1=1.19 

80  SUM=SUM*H!I) 

AVE  = SUM/  19 

CALL  aALHGT!H9.N9.VD. 1.15.1 .2.. 01 . AVE . CM . TM . H . BAL ) 
85  CONTINUE 

IF  ! IPR1NT.lt. 1 ) GO  TO  86 
PRINT  9000 
PRINT  9070 
PRINT  9062 

PRINT  901 0. !H! I ) . f =1 . T 9) 


PRINT  9064 

ORINT  9010.( T(  I ) « 1 = 1 « I 9) 

PRINT  9066 

DRINT  901 0, <U( I ) . I =1 . 19) 

PRINT  906fl 

PRINT  9010. (V(  I ).  1 = 1 * 19) 

^16  CONTINUF 

IF  (IFOR.UT.l)  GO  TO  95 
C 

C UPDATE  ALL  GWC  FORECAST  FIELDS  BY  AOVFCTING  WITH  COMPUTED  TRENDS 

r 

DO  90  IH=IG.i;»0 
N 1=  ( IM-  1 ) *10*(  L-  1 ) *K  3AN494-  1 
N2=Nl ♦NAR-l 
DT=DTS 

LIH=1  lH-1  )*L4FL 
IF  (IH.GT.IG)  GO  TO  SB 
DT=-FG 1 * I 0. 0*OT 
1 = 1 

DO  87  N=N1 ,N2 

U( I )=OATA(N )*FG1+0ATA( N+ID) *FG2-U( I ) 

V(I)=0ATA(N  + N4  9)«‘FG1*DATA(N^N4  9-HD)*FG2-V1I) 
IHl=lNT(DATA<N+N98  >*0.0001) 

IH2=INT(DATA(N+N98-*-ID  ) *0.0001  ) 

I Tl=OATA(N*N98  )-IH!*10000 
IT  2=0  AT  A { N + N98-HD  )-IH?*10000 
I T1=IA8S{ ITl ) 

IT2=IABS( IT2) 

H(I)=  IHl  *FG1*IH2*FG2-H(  I ) 

T(I)=  I T 1*FG 1+ 1 T2*FG2-T( I ) 

87  1 = 14-1 

88  LIH=(  IH-1  )*L44-L 

CALL  ADVECI DT.TRU(LIH).TRV(LIH) .U.BAL.CM) 

CALL  AOVECl OT. TRUIL IH ) .TRV(LIH).V. BAL.CM ) 

CALL  ADVECIDT.TRUILIH)  . TR V ( L 1 H ) .H . HAL . CM ) 

CALL  AOVFCl DT. TRUIL IH) .TRV(L I H ) . T . B AL , C M ) 

1=1 

DO  89  N=NI.N2 

OATAl  N)=OATA(  N)-U(I)*10 
OATAlN+NARlsDATAl  N4-N49)-V  ( I ) *1  0 
IHl  = INT(OATA(N4-N98  )*0.0001  ) 

I T1=DATA  ( N4-N98  )-IHI*10000 

IT  1= lABSI  ITl  ) 

IH2=IH1-H(1)*I0.0 
IT2=IT1-T(l)*l0.0 
IF<  IH2.lt, 0)  IT2=-IT2 
DATA!  N4-N98)=IH2*1  00004- IT2 


89  1=1*1 

90  CONTINUE 
95  CONTINUE 

r 

C FNO  OF  LEVEL  LOOP 

C 

r 


BEST  AVAILABLE  CCPY 


~X- 


n r, 


C WRITF  OUT  UPDATAO  FIFLOS  ON  TAPES 
C 


IF  (ITAPE.LT.l)  GH  TO  100 
WRITE  (3)  JDATE.JTIMF 
DO  9fl  1H=1 • 120 
DO  9S  L-Lt«l.4 

Nl  = < IH-1  »*IO<-IL-l  )*K3*N49-H 
N2=Ni  + N49-l 


98  WRITE  (3)  (OATA(N)  .OATA<N  + N49)*r)ATA(N>N98)  tN=Nl  .N2> 

100  CONTINUE  ) 

END  OF  ARMY  LOOP  1 

101  CONTINUE  ^ 

GO  TO  1000  jj 

999  PRINT  9999  ii 

, lOOO  CONTINUE  I 

9000  FORMAT!  1 HI  » > 

[ 9001  FORMAT! 2X ,8HCARO  A ,TI10//>  f 

9002  FORMAT!  2X,aHCAR0  B . 2 I 1 0 , •SF  1 O . 2// > 

9003  FORMAT!  2X»  8HCARD  C .110//)  ij 

9004  F0RMATI8X, 5110/1  ) 

9005  FORMAT  l/2X,eHCARO  E tllO*//) 

9005  FORMAT! 10X«I5»SX.3F10.?/»  I 

9007  F0RMAT!1AX» I2«2X.4Fi0.2)  3 

9010  FORMAT!  10X«7F5.0)  Sj 

9025  FORMAT!  IH  ) ji 

9050  F0RMAT!30X,27HINPUT  DATA  FOR  PDRR  ROUTINE///)  1 

9054  FORMAT!  2X.5THCAR0S  O YEAR  MONTH  DAY  HOUR 

2 MlN/l  I 

9055  F0RMAT!2X,8HCARD5  F , 9X . 20 HS T AT  I ON  I NFORMAT I ON// 1 2X , 1 0 HNUMBER  j 

2 .30HL0NGITUTE  LATITUTE  FLFVATION/)  S 

9057  FORMAT!25X. 13HSTATION  DATA//  9X.33H  NUMBER  O VALUF  TEMPERAT  j 

2URE.  19H  U COMP  V COMP  ./)  j 

9059  FORMAT! 10Xt22HOE5ULTS  FOR  LEVEL  L * .11//)  ’ 

9050  FORMAT! 20X.21HGWC  DATA  ON  MINI  GRID//)  j 

9052  FORMAT I/24X.8HD  VALUES/)  ! 

9054  F0RMATI/23X.11MTEMPERATURE/)  ' 

9055  FORMAT! /25X.6HU  COMP/)  ? 

9050  FnRMAT!/25X »5MV  COMP/)  j 

9070  FORMAT! I 5X. 29HUPDATEO  GWC  DATA  ON  MINI  GRID//) 

9100  FORMAT I///45X.34HPR0GN0ST IC  DATA  REANALYSIS  ROUTINE.////// 

2///.  57X.5HUN1 TS./// .53X .1  IHSPEFD  - MRS . // , 53X . 1 5HD I RFCT I ON  - OE  i 

3G,//.  53X. 15HHF IGHT  - ME TF RS. // . 53X . 1 9HTEMOE RA TURE  - OFG  K» 

9999  FORMAT! 10X.23MINCONSI5TENT  INPUT  DATA) 

STOP 
END 
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SUBROUTINE  GtVCINd  . 1 G WC  • JOATE  . JT  I ME*  OAT  t • IMOI 
C 

C THIS  SUBROUTINE  READS  IN  THE  GMC  FORECAST  DATA. 

C 

C 1 = HOUR  INDEX 1 EQUALS  1 FOR  INITIAL  ANALYSIS  TIME  AND  I EQUALS  5 

C FOR  FORECAST  4 HOURS  AFTER  INITIAL  ANALYSIS 

C IGWC  = AO  TIME  IN  MINUTES COMPUTED  IN  SUBROUTINE 

C OATl  = ARRAY  CCNTAINING  GWC  DATA 

C IMD  = NUMBER  OF  DAYS  FROM  BEGINING  OF  YEAR  UP  TO  MONTH  INDICATED 
C BY  INDEX 

C 

DIMENSION  I ATI  600)  • IPREI 5 ) «PRS< 5) 

DIMENSION  OATI  I 1 )•  IMD(  1 ) 

DATA  RH  /50.0/ 

DATA  LI .LA.KA.NZS.LlOO/l.A.A.ZS* 100/ 

DATA  (PRS=870. 0*850.0. 700.0.500.0. 300.0) 

IF  II.GT.l)  GO  TO  to 

READ  (I)  NO. (IPRE(N)*N3l*S) 

PRINT  I .1PRE12) *1PRE(3 ) .IPREIA) .IPREI 1 ) 

I FORMAT! 10X.42HGWC  FORECAST  DATA  READ  IN  FOR  ANALYSIS  OF  *I5*IH/*I2 
2 .1H/.I2.5X. IS.4H  GMT//) 

IPREtsINTI IPREI l)*O.Ot ) 
lPRE3=IPREt3) 

IGWC=IPREI1  )-I PRE 14 40 +1440*1 IPREI 2) ♦IMDI  I PRE3 ) +365* IPRE I 4 )- 366 ) 
JOATE=  IPRE  14  )*I0000+  IPREI  3 )•  1 00  + 1 PREI  2) 

JTIME=1PREII ) 

10  READ  I I ) lAT 

IN=0 

DO  20  L=L1*L4 
LT=IL- I )*N25 
DO  20  K=1.K4 
LK=IK-1)*L100+LT 
DO  15  N=1*N25 
LN=LK+N 
IN= IN+l 
XAT=I ATILN) 

IF  (K.NE.K4)  GO  TO  15 
XAT=XAT-273.16 

yat=xat 

RHT=RH*0.0I 

IF  IL.EQ.l)  RHT=1 ATILN+LIOO )*0.0I 
EXN=7.5*XAT/I XAT+237.3) 

EST=6*1 1*10.0 **EXN 
RMaRHT*0.622*EST/IPRS<L)-EST ) 

XAT= IX AT +273.16)*! 1 .0+0.6 1*RM) 

IF  IN.NE.13)  GO  TO  15 
15  OATH  IN)=XAT 
20  CONTINUE 
RETURN 
END 


46 


•X' 


SUBROUTINE  L VCONVI U • V .OLONI 


C THIS  SUBROUTINE  CONVERTS  CWC  GRID  U ANO  V COMPONENTS  INTO  SPHERICAL 
C GRID  U ANO  V COMPONENTS. 

C 

C U.V  = INITIAL  VALUES  ARE  C«C  COMPONENTS 
C = FINAL  VALUES  ARE  SPHERICAL  COMPONENTS 

C OLON  s LONGITUDE  WEST  OF  GREENWICH 
C 

SPO*SQRT( UWU*V*V) 

ALPH=  ATANtU/(ABS(V)«^0.0000  1 ) >*5.72968 
IF  (VI  1.2.2 

1 ALPH=36.0+ALPH 
GO  TO  3 

2 ALPH=18.0-ALPH 

3 ALPH=ALPH-0L0N*0 .1 *8.0 
OIRsALPH*t8.0 

IF  (OIR.GT.36.0)  OlR=OIR-36.0 
ANG«OIR*0. I 74S33 
U=-SPD*S INI ANG  > 

V=-SPD*COSI *NG> 

RETURN 

ENO 


) 


r>  r\ 


SUUROUT I NF  |NTPT( I . 1 R»  A* OX . OV • AS I 


THIS  SUBROUTINE  INTERPOLATES  A VALUE  UETNEEN  FOUR  GRID  POINTS. 

C 

C I ■=  UPPER  left  grid  POINT  OF  THE  FOUR 
C IR  s NUMOER  OF  POINTS  IN  A RO« 

C A = ARRAY  CONTAINING  THE  GRID  POINT  VALUES 

C OX.OY  * FRACTICNAL  distance  BETVEEN  grid  points  AT  AHICM  INTERPQLATfO 
C VALUE  IS  TO  BE  OBTAINED 

C AS  « INTERPOLATED  VALUE 
C 

DIMENSION  AC  II 
X1  = AI  I ) 

X23 A( I ♦ I I 
X3*AI  I^IR  I 
XAsAl  I^IM'fl  I 
A I«X1 
A3>X2 
A7«X3 
A9SXA 

A2stX2>XI  ) AOXFXI 
AA^I  X3-XI  l*OV>  XI 
A6z( XA-xalAOV^X? 

A 8=  I XA-X3)*OX«'X3 

AS>  I OX  *4  AA-A4  t ♦A«F0TA(A8-A2>«’A2)A0.S 

RETURN 

END 
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r>n  nr>  n nr>r>n  n nn  or^r^nr^nnn^nn^>r^^inr> 


SUBROUTINE  T RENO ( MVARS «NS I 2E « OX . DY . OAT  A t . DAT A2.  16) 

C CALCULATE  THE  DISPLACEMENT  NECESSIARY  TO  GIVE  THE  SMALLEST  ROOT  MEAN 

C SQUARE  DIFFERENCE  BETWEEN  CORRESPONDING  DATA  VALUES  OF  OATAI  AND  DATA2 

C 

C FORMAL  PARAMETERS  ............ 

C MVARS  NUMBER  OF  VARIABLES  (MAXIMUM  IS  10) 

NSI2E  GRID  SIZE  ASSUMED  TO  BE  THE  SAME  FOR  X AND  V DIRECTIONS 

DX.OV  ON  INPUT.  THE  INITIAL  GUESS  FOR  THE  TREND  IN  THE  AREA*  ON 
OUTPUT.  THE  CCMPUTEO  TREND  (TREND  IS  GRID  DISPLACEMENT 
FOR  BOTH  CASES). 

OATAI  AND  0ATA2  ARE  THE  ARRAYS  OF  DATA  AT  TIME  I AND  TIME  2 (OR 

LOCATIONS  I AND  2).  RESPECTIVELY.  BOTH  ARE  ARRAYS.  NSI7E 
BY  NSI2E  (GRID  SIZE)  WITH  NVARS  VARIABLES. 

THESE  DATA  ARRAYS  ARC  ASSUMED  TO  BE  THREE  DIMENSIONAL. 

PACKED  BY  (X. Y.VAR  lABLE  ). 

LOCAL  PARAMETERS  ................. 

NTEMF  TEMPLATE  SIZE  (IN  GRID  UNITS).  BOTH  X AND  Y DIRECTIONS 

PTSMIN  MINIMUM  NUMBER  OF  POINTS  FOR  MAKING  A CALCULATION 

OXYMAX  INITIAL  GUESS  MOTION  — (OX.OY)  ON  INPUT  — MUST  BE  LESS  THAN  THI 
MAXVAR  MAXIMUM  NUMBER  OF  VARIABLES 


logical  debug 

DIMENSION  OATAI ( NS I ZE • NS  1 2E .MV ARS I * O AT A2 ( NS  I ZE. NS I ZE .M V ARS ) 

01  ME NS  I ON  RMS ARY (5.5) * AVEl ( 10 ).S0V1( 1 0 ) . AVC2(  10) . S0V2(  10) 

DIMENSION  IL0C(9). JL0CI9) 

INITIALIZE  PROGRAM  CONSTANTS 

DATA  DEBUG/. FALSE./*  LP/6/*  MAXVAR/IO/.  OXYMAX/5.0/*  PTSMIN/IO.I/ 
ARRAYS  FOR  PICKING  THE  NEXT  LOCATION  TO  MAKE  A CALCULATION 
DATA  (LOC/O. -1.0.1 .0.-1 *1 *1 *-!/■  JLOC/ 0 . 0 . 1 . 0. - 1 . I . 1 * - 1 * - 1 / 


INTERPOLATE  THE  LOCATION  OF  THE  MINIMUM  RMS  DIFFERENCE.  GIVEN  THE 
APPARF^'  MINIMUM  AND  THE  VALUE  ON  EACH  SIDE  OF  THE  MINIMUM 
(FINOS  . <£  LOCATION  OF  THE  MINIMUM  OF  A PARABOLA  THRU  THESE  3 POINTS) 
OMINlRI.R2.R3)  = O.S*(R3  - R1 ) / ( R2 *2 . 0-R1-R3 ) 

the  minimum  RMS  DIFFERENCE  AT  THIS  LOCATION  IS  GIVEN  BY  THIS  EQUATION 
RMINIR1.R2.R3)  = R2  * ( R3-R 1 ) A*2/ ( B . 0* I R2*2 .0 -R3-R I ) ) 


INITIALIZE  SUBROUTINE  PARAMETERS 
NTEMP  s MINO(NSIZE.IO) 

NVARS  > MINOIMVARS. MAXVAR  ) 

YF  40EBUG)  PRINT  1001.  NS  I ZE* NV ARS. DX. OV 

MAKE  S-'Rc  INITIAL  ESTIMATE  OF  MOTION  IS  WITHIN  BOUNDS 

IF  fA8S(0X).LT. OXYMAX  .ANO.  ABS( OY ) .L T. OX YM AX ) GO  TO  110 


XR  - 

OX 

VR  = 

DY 

DY  = 

AMAXK-OXYMAX* 

AMINIIOY. 

OXYMAX) ) 

OX  s 

AMAXl (-OXYMAX. 

AMINl (OX* 

OXYMAX ) ) 

PRINT 

1000*  XR.YR 

•OX.OY 

110  CONTINUE 

C CALCULATE  MEAN  ANO  STANDARD  DEVIATION  OF  EACH  VARIABLE 

C THESE  ARE  USED  LATER  IN  THE  PROGRAM  TO  SCALE  THE  DATA  FOR  EACH 

C VARIABLE  SO  DIFFERENCES  OF  ALL  VARIABLES  (tONTRIBUTE  EQUALLY  TO  THE 

— ■ ■ - ■ ■ 

This  subroutine  was  written  by  Dan  Wolf  (SRI) 
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U U 


ROOT 


MEAN  SQUARE  DIFFERENCE. 

DO  140  LV=I«NVARS 
AVEt(LV>  = 0.0 

AVE2<LV}  = 0.0 

SOV2(LV)  = 0.0 

SOVULV)  = 0.0 

OO  130  Isl.NSIZE 
DO  130  Jsl.NSIZE 


AVEl (LV) 

= AVEIILV) 

♦ 

OATAl ( 

AVE2(LV) 

= AVE2(LV) 

♦ 

DATA2( 

SOVl (LV) 

= SDVKLV) 

♦ 

OATAl ( 

S0V2(LV) 

= S0V2(LV) 

4- 

DATA2  ( 

I • J.LV) **2 


I 30 
1 40 


continue 

CONTINUE 

PTS  = FLOAT(NSIZE4NSIZE> 
DO  150  LV=1«NVARS 


1 50 


AVEl(LV)  = 
AVE2CLV1  = 
SDVKLVI  = 
SDV2(LVI  = 

IF  <SOVt(LV) 

IF  (SDV2<LV) 
CONTINUE 

IF  (DEBUG!  PRINT 


AVEl (LVl/PTS 
AVE2(LV1  ✓PTS 
SORT ( (SOV  1 (LV ) 
SQRTI  ( SOV2(LV) 
.LT.  0.011 
.UT.  0.01) 


- PTS*AVE1 (LV)**2) 

- PTSAAVE2(LV)**2) 
SDVKLV)  = 1.0 
SDV2CLV)  = 1.0 


CPTS-  1. 0)  » 
(PTS-1 .0) » 


152  . 


152  FORMAT  I44H0VAR  AVEl 

* •/.( 1X.I3,4F10.2) ) 

INITIALIZE  ARRAY  OF  RMS  DIFFERENCES 
OO  170  1=1.5 

OO  170  J=1.5 

RMSARVd  .J)  = -1.0 

170  CONTINUE 


(LV. AVEl (LV) .SOVl (LV) 

.AVE2(LV) .SOV2(LV)  .LV=1  .NVARS) 
SDVl  AVE2  SOV2 


OF  THE  AREA 


C INITIAL  SECTION  OF  *DATA1*  FROM  CENTER 

LIl  = ( NSIZE-NTEMP)/2  ♦ 1 

LJl  = ( NSlZE-NTEMP)/2  ♦ I 

C DISPLACE  THE  SECTION  FRCM  40ATA2*  BY  THE 
C MOTION  FOR  THE  AREA  (OX.OY) 

C THIS  IS  THE  LOCATION  OF  THE  INITIAL  MINIMUM 
' IMIN  = MAXO (LI !♦ IFIXIOX ) . 1) 

JMIN  = MAX0(LJ14IFIX(DY).  1) 

C INITIALIZE  THESE  PARAMETERS  EACH  TIME  THIS  SUBROUTINE 

NLOC  = 0 

BEST  = 99999.9 

LMIN  = 9 


INITIAL  GUESS  OF  THE 


IS  CALLED 


LOOP  TO  TRY  AT  LEAST  3 SUCCESSIVE  LOCATIONS  IN  BOTH  X AND  Y DIRS. 
200  NLOC  = NLOC  * 1 

IF  (NLOC  .GT.  9 ) NLOC  = 2 

IF  (NLOC  .EQ.  LMIN)  GO  TO  300 
210  L12  = 1M1N4^  1 LOCI  NLOC) 

LJ2  = JMIN  4 JLCC(NLOC) 


C 

IF  (DEBUG)  PRINT  221.  NLOC. L I 1 .L J 1 . L 1 2.L J2. NS IZE. NTEMP 

221  FORMAT  I 9HOLOCAT I ON . 1 2 . 9H  IN  DA TAl . 2 IS .SX . 8H I N OATA2.2I5.5X 

♦ .9HGRIO  SIZE. I5.5X. 13HTEMPLATE  SIZE. IS) 

C SEE  IF  CALCULATION  HAS  BEEN  MADE  AT  THIS  LOCATION 
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s LI2  - IMIK  ♦ 
* LJ2  - JMIN  ♦ 
(RHS«RY( II.Jl I .GT. 


GO  TO  200 


COMPUTE  THE  CORRELATION  BETWEEN  THE  NTEMP  BY  NTEMP  DATA  POINTS. 
ITOP  LEFT  CORNER)  AT  (LIl.LJl)  OF  DATA  1 AND  (L12.LJ2)  OF  i>ATA2. 
GET  BEGINNING  AND  ENDING  LOCATIONS  IN  X AND  Y DIRECTIONS 


IE  = MINOINTEMP.  NS  I ZE4- 1-L 1 1 . NS1ZE4-1 

IB  = MAXOI  1 . 2 -Lit*  2 

JB  = HAXOI  1 • 2 -LJl.  2 

JE  = MINOINTEMP.  NSIZE«1-LJ1.  NSIZE^l 

PTS  = FLOAT(IJE-JB*l )•< IE-IB+1 )*NVARS) 

IF  (PTS  .LT.  PTSMIN)  GO  TO  200 
FORM  SUM  OF  SQUARED  DIFFERENCES 
RMS  = 0.0 

00  250  LV=1.NVARS 

SDV  = 0.5  * (SOVl  (LV)4-SDV2(LV)  ) 

OO  240  I=1B. IE 


NSIZE4  1-LII. 
2 -L  I 1 , 
2 -LJI. 
NSIZE41-LJ1. 


NSIZE  41-^2  ) 
2 -LI2) 

2 -LJ2) 

NS IZE41-L J2) 


240 

I=IB.  IE 

I 1 = 

I ♦ LIl  - I 

12  = 

I ♦ LIZ  - 1 

00  230 

j=ja. JE 

J1 

= J ♦ LJl  - 

1 

J2 

= J ♦ LJ2  - 

1 

RMS 

= RMS  * ((OATAKIl 

RMS  = RMS  ♦ I (OATAl < I 1 , J1 ,LV )-0ATA2 I 12. J2.L V) )/SDV )4*2 
) CONTINUE 

} CONTINUE 

3 CONTINUE 

CALCULATE  ROOT  MEAN  SQUARE  DIFFERENCE 
RMS  = SQRTIRMS  / PTS) 

II  = LI2  - IMIN  ♦ 3 

J1  = LJ2  - JMIN  ♦ 3 

RMSARYIIl. Jt ) = RMS 

IF  (RMS  .GT.  BEST)  GO  TO  290 

NEW  MINIMUM  RMS  VALUE  HAS  BEEN  FOUND  MOVE  POINTERS  TO  NEW  LOCATION 

IF  (NLOC  .EQ.  1)  GO  TO  280 

SHIFT  4RMSARY*  VALUES  SO  THE  MINIMUM  VALUE  IS  AT  (3.3) 


lOIF 

= 

1 LOCI NLOC) 

I INC 

1 ♦ 

IDIF4< 

1 

- 

IDIF) 

IB 

= 

1 * 

I0IF4(- 

2 

♦ 

24IDIF) 

IE 

5 ♦ 

IOIF*( 

1 

- 

2*IDIF) 

JOIF 

= 

JLOCINLOC  ) 

JINC 

1 ♦ 

JOIF*( 

1 

- 

JDIF) 

JB 

1 ♦ 

JDIF*|- 

2 

♦ 

24J0IF) 

JE 

= 

5 ♦ 

J0IF4< 

1 

- 

24JDIF) 

I 

IB  - 

I INC 

J 

= 

JB  - 

J INC 

J 

s 

J 

* JINC 

J1 

J 

4 JOIF 

I 

= 

I 

♦ lINC 

11 

s 

1 

♦ IDIF 

RMSARY(I.J)  s 
IF  (I  .NE.  IE) 

= IB  - I INC 


RMSARYI II.  J1  ) 
GO  TO  262 


n n nnnn  on 


c 


IF  (J  .NE 


JE)  C-O  TO  260 


263 

264 


265 

266 


280 


281 


IF  (JOIF  .EQ.  0)  GO  TO  264 
J = JE  ♦ JINC 
DO  263  1=1.5 

RMSAMYll.JI  = -I.O 
CONTINUE 

IF  (lUIF  .Ed.  0)  GO  TO  266 
I = IE  ♦ IINC 
DO  265  J=l.5 

RMSARV(I.J)  = -1.0 

CONTINUE 
CONTI  NUE 

SAVE  THE  MINIMUM  VALUE  AND  ITS  LOCATION  tlN  OATA2) 

LMIN  = NLOC 
BEST  = RMS 
IMIN  = LI2 
JMIN  = LJ2 

IF  lOEBUG)  PRINT  281.  I < RMSAR Yl I . J I . I = 1 ,5  I . J= 1 « S t 

FORMAT  (10X.5FI0.3> 


READY  TO  PROCEED 
290  CONTINUE 

IF  (DEBUG) 


♦ 


TO  THE  NEXT  LOCATION 

PRINT  291.  NLOC. LMIN. IMIN. JMIN 

. IB. IE. JB. JE.PTS.  RMS 


291  FORMAT  ( 1 OX. 81 5 .F6 . 0 .F 1 5 .6 ) 

IF  (LMIN  - NLOC)  200.210.200 


MAKE  SURE  THAT  A CALCULATION  HAS  BEEN  MADE  AT  LOCATION  9 OF  SEARCH 
300  IF  (NLOC  .EQ.  9}  LMIN  = 2 

IF  (NLOC  .EQ.  9)  GO  TO  210 

INTERPOLATE  EETWEEN  THE  3 VALUES  IN  BOTH  DIRECTIONS  TO  GET 
LOCATION  OF  MINIMUM  RMS  DIFFERENCE 
OX  = 0.0 

OY  = 0-0 

IF  (RMSARY(2.3) .LT.0.0  .OR.  RMSAR Y< 4 • 3 ) .LT . 0 • 0 ) GO  TO  310 
OX  = 0M1N(RMSARY( 2. 3 ).RMSARY( 3. 3) .RMSARYI  4. 3) ) 

XR  = RMIN(RMSARY(2.3I .RMS  ARY ( 3.3>.RMSARY(4.  3) ) 

310  IF  (RMSARYI 3. 2>.LT.O. 0 .OR.  RMSAR Y( 3. 4 ) .L T. 0. 0 ) GO  TO  320 

OY  = DM1N(RMSARY(3.2),RMSARY(3.31.RMSARY(3.4)) 

YR  = RM>IN(RMSARY(  3. 2>.RMSARY(  3.3)  .RMSARYI  3.4)  ) 

320  DX  = OX  ♦ FLOAT!  IMIN-L It  ) 

DY  =-(OY  ♦ FLOAT! JMIN-LJl )) 


PRINT  IF  DEBUGGING  FLAG  IS  SET 
IF  I, NOT. DEBUG)  RETURN 

PRINT  1002.  ( I.  I=I.S). ( J. (RMSARYI I . J) . 1 = 1.5) . J=I . 5) 

PRINT  1003.  OX.OY.XR.YR 

RETURN 

1000  FORMAT  (I9H01N1TIAL  ESTIMATE  ( .F5. 2 . t H. .F5. 2 

*■  .29H)  OF  (OX.OY)  IS  TOO  LARGE.  ( .F5. 2.  IH.  .FS.2 

♦ .IIH)  HAS  USED.  ) 

1001  FORMAT  ( I7HOSUBROUTINE  TREND. 2 1 5. 2F 1 0 .3 ) 

1002  FORMAT  (22H0TREN0  RMS  0 1 FFERENCES . / « SX. 5 1 10 . / . ( IX • I 4. SF 1 0. 5) I 

1003  F0RMAT(9X. 12HOISPLACEMENT.2F8.3.9X.13HMIN  RMS  0 I FFS « 2F 1 0 . 6 ) 

END 
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SUBROUT INE  RA«1N(  JT. I .ES • UR • VR • HR  « TR « I MO. lOATE . IDIF. I GMT) 

C 

C THIS  SUBROUTINE  READS  IN  THE  U S ARMY  RAWINSONOE  DATA. 

C 

DIMENSION  UR( 1 ) . VR( 1 ) .HR!  1 ) .TR(1 ).ES( 1 ) . IMO< 1 ) 

DIMENSION  OAT(5t2)«ISi25).HOZ(  17) .PR SI  6) .STOf 6) 

DATA  ACR.CKM.CC.XNIE/O  .01  74533*  0.  5 1 48 . 29 . 29.-<)99. 9/ 

DATA  J10.J17.J24.J28.JS1.K4/10.17.24.28.5I.4/ 

DATA  STD/  0.0. 1457. 0.3012.0.5574. 0.0. 0.0.0/ 

DATA  PRS/  870.0.850.0.700.0.500.0.300.0.200.0/ 

DATA  HQZ/0.0. 100 . 0 . 350. 0 • 750. 0 . 1 250.0. 1 750.0.2250.0.2750 .0 .3250 .0 
2 .3750. 0.4250.0.4750.0.5500.0.6500.0. 7500.0.8500. 0.9500.0/ 

IMIN=1440 
ESST=l 300.0 
DO  18  J=I.JT 

18  1S(J)=0 
JK=  JT*K4 

DO  19  J=1«JK 
HRIJ )=XNIL 
TR(  J)  = XNIl. 

UR( J)=XNIL 

19  VRCJ)=XNI1. 

IF  I I.GT.l  ) GO  TO  22 
C 

C READS  IN  DATA  SET 
C 

20  READ  2.0AT 

2 FORMAT!  10F8.  1) 

C 

C LOOP  THROUGH  DATA  SET  IN  RECORD 
C 

DO  30  J=1.J10 
22  JC=I J-1 )* J5I+1 
C 

C CHECKS  DATA  FOR  DATE . T I ME « AND  TYPE 
C 

lPl=INT(OATf JC)41.0E-4) 

DATJ=0AT(JC)-IP1*1 .0E4 
IP2=INT(DATJ*l .Oe-2 ) 

IP3  =DATJ-IP241.0E2 
IP4-INT<0AT( JC41 )*0.01  ) 

1P5=0AT<  JC+I)-IP4*100 

JDATEs IP5460*( IP4424*< lP34lM0f 1 P2 ) +3654I P 1-366) )-IGMT 
IF  ( JDATE.LT.IDATE-IMIK)  GO  TO  32 
IF  ( JDATE.LT.IOATE-IDIF)  GO  TO  30 
IF  I JOATE.GT.IOATE4IOIF)  GO  TO  32 
JS=INT(DAT( JC42)*0. 1 ) 

IT  =0AT<JC42)-JS410 
IF  ( 1T.LT.4.0R.IT.GT.6)  GO  TO  30 
IF  < IS(JS).EQ.6)  GO  TO  30 
NOATE-lNT(OAT<  JC) ) 

NTIME=INT<OATI JC+1 )) 

ISI JSI^IT 
ES3S=ES(  JS) 

C 
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on  on 


SECECTS  OUT  HEIGHTS  AND  TEMPERATURES  FOR  THE  MANDATORY  UEVELS 


K=l 

DO  26  JQ=4,J24.4 
JK=(K-1I*JT  ♦JS 
JCQ-JC4 JQ 

IF  (K.EQ.t)  GO  TO  25 

IF  <K.GT.K4)  GO  TO  27 

IF  (DAT! JCQ) .NE.PRS(K) ) GO  TO  26 

25  IF  {OAT ( JCO-1 I .LT .0 .0 > GO  TO  125 
HR«  JK)=UAT{ JCQ-I l-STO( Kl 

TR( JKI^DATI JC041 >4273. 16 
IF  (K.NE.I ) GO  TO  125 

PRSI=I lESJS-ESST) /«TR( JK) ♦CC)4I .0  »*DAT<  JCQ» 

HR (JK )«-( DAT! JCQ  )-  PRSl  » ATR ( JK I ACC /OAT( JCQ ) 

125  K=K4| 

26  CONTINUE 

27  KT=K-l 

INTERPOLATES  ZONAL  WIND  VALUES  TO  OBTAIN  VALUES  FOR  MANDATORY  LEVELS 
JC J=JC4J28-1 

IF  ( ABS (DAT IJCJ )) .GT .360.0 ) GO  TO  127 
ANG=ACR«OAT( JCJ) 

UR( JSI=-OAT( JCJ41 ) ♦SINI ANGI ACKM 
VRI JS»=-DAT ( JCJ41 ) *COS( ANG) 4CKM 

127  J2=2 

DO  29  K=2.KT 
JK=(K-n*JT  4JS 
JG= JCJ 

DO  28  JHSJ2.J17 
JG= JG42 
JI  = JG4-I 

HRJK=HR(JK)4ST0(K) 

IF  ( HOZ(JH)4ESJS.LT.HR JK)  GO  TO  26 

RH=  (HRJK-HOZ( JH-l ) -ES JS ) / I HOZ IJ HI-HOZ ( JH- I I ) 

IF  (ABS(OAT(  JG) I .GT.360. 0)  GO  TO  128 

UR( JK)  = <-RH*OAT ( J I ) *S I N ( ACR *0 AT ( J G ) >40 ATI J 1 -2 ) *SI N( ACR4 DA T( JG-2 I ) • 
2(RH-1. 0) I «CKM 

VM(  JK»=(-RH*DATC  JI  )«<COS(  ACR*DAT(  JG)  ) 4D AT  ( J I -2 ) ♦ COS(  ACR*OAT(  JG-2>  )♦ 
2 (RH-1.0 ) I ♦CKM 

128  J2=JH 

GO  TO  29 

28  CONTINUE 

29  CONTINUE 

30  CONTINUE 
GO  TO  20 

32  CONTINUE 
JSM=0 

DO  150  J^l.JT 
IF  (ISIJI.GT.O)  JSM=JSM4l 
150  CONTINUE 

PRINT  l.JSM 

1 FORMAT! 10X.37HU  S ARMY  RAWINSONOE  DATA  READ  IN  FOR  • 1 2 1 9H-S TAT I ONS 
2 .//) 

RETURN 

END 
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SUBROUTI  Nf!  MESH(  KQ.  J J.  YS«  XS  t VS.US>HS*  TS  .YLtXL.VL  .UL.HLt  TL  I 


C 

C THIS  SUBROUTINE  COMPUTES  GRID  POINT  VALUES  FROM  OBSERVED  DATA  AND  AN  INITIAL 
C GUESS  GRID  POINT  FIELD  I IF  KQ  = 2)  BY  A LEAST  SQUARES  FITTING  OF  THE  DATA 

C (SEE  ENDLICH  AND  MANCUSO.  MON  NEA  REV.  1968*  342-3501.  IT  CAN  ALSO  BE  USED  TO 

C GENERATE  ITS  DUN  INITIAL  GUESS  FIELD  IIF  KQ  = I).  ALTHOUGH  THIS  PARTICULAR 
C SUBROUTINE  HAS  BEEN  SET  UP  FOR  ANALYZING  MIND  COMPONENTS.  IT  CAN  ALSO 
C BE  USED  TO  ANALYZE  ANY  3 QUANTITIES  THAT  ARE  PLACED  IN  ARRAYS  HS.US. 

C AND  VS  (IDS  SHOULD  BE  SET  TO  ZERO  IN  ORDER  TO  BYPASS  THE  UPSTREAM. 

C TIME.  AND  HEIGHT  WEIGHTING  EFFECTS) 

C 

C THIS  SUBROUTINE  ALSO  INCLUDES  THE  ENTRY  POINT  TO  SUBROUTINE  MESHSET. 

C 

C JJ  = NUMBER  OF  WIND  DATA 

C KS  = NUMBER  OF  CLOSEST  DATA  TO  A GRID  POINT  USED  TO  COMPUTE  ITS  VALUE 

C Wl  = WEIGHT  GIVEN  TO  INITIAL  GUESS  VALUE  AT  GRID  POINT 

C C2  = WEIGHTING  CONSTANT 

C APHl.APH2.APH3  = WEIGHTING  FACTORS  FOR  UPSTREAM,  TIME.  AND  HEIGHT  EFFECTS 
c RMAX  = MAXIMUM  RADIUS  OVER  WHICH  THE  KS  CLOSEST  DATA  ARE  SELECTED  FOR 
C ' COMPUTING  A GRID  POINT  VALUE  lOEG  LATITUDE) 

C X15  = LONGITUDINAL  RANGE  I+XIS  TO  -X15)  OVER  WHICH  THE  SEARCH  FOR  THE 
C KS  CLOSEST  DATA  IS  MADE  (DEG  LONGITUDE) 

C YS.XS  = LATITUDE  AND  LONGITUDE  OF  WIND  DATA  IDEG) 

C US* VS  = U AND  V COMPONENTS  OF  WIND  DATA  (M  SEC-1) 

C UN.VNIKQ=1)  = INITIAL  VALUES  NOT  USED 

C = FINAL  VALUES  A SMOOTH  ANALYSIS  USED  AS  THE  INITIAL  GUESS 

C field  IN  THE  KQ=2  COMPUTATION 

C UN.VN(KQ=2)  - INITIAL  VALUES  USED  AS  THE  INITIAL  GUESS  FIELD 

= FINAL  VALUES USED  AS  THE  FINAL  ANALYSIS 

HS.HN  = INITIAL  DATA  AND  ANALYZED  GRID  POINT  VALUES  FOR  AN  ARBITRARY  QUANTITY 
C 

DIMENSION  OSI  SO > * JS I SOI.IOI  S0).OVR(20) 

DIMENSION  YSI  1 ) .XSI  1 ) . VS(  1 ) *US( I ) .HSI I ) .TS( I ) 

01  MENS  ION  YLI  1 )*  XL(  1 )•  VLI  1 I.ULI  I I.HLI  1 ) . TLf  1 ) 

COMMON/CIS/  ISIIOOO) 

COMMON /C  GO/  M9.N9.I9*Me*N8.YB.XB.OO 
COMMON/CPTS/  KS.WI .C2.RMAX.KSS5. lOS.KSW.ALPH 

EQUIVALENCE  lOVRf 1 ) . ONH ) * I O VR ( 2 ) . DHH ) . (OVR(3) .OUH) . I DVR I 4 ) • OVH ) . 

2 <0VR(5).0TH) . (DVR (6) .OXH). (OVR( 7 )*OYH). 1 0 VR I 8 ) . OXTH ) . 

3 I0VR(9).0XXH) .IDVRI 10) .OVVH) . lOVRI 1 1 ) .OXHH).  (DVRl 12).  OXUH). 

4 IDVRI  13)*OXVH).  (OVRI  14)  .OYHH)  * (DVR  I 1 S)  .OTUH  I . IDVRI  16)  .OWN)  . 

5 IDVRI 17 ) .OXTH) . (OVRI 18) .OYTH) 

IF  (KQ.LT.l)  GO  TO  18 
K0S=K0-14KSS5 

MsO 
L = 0 
1 = 0 

75  M=M+l 

IF  (M-M9)  77.77.100 
77  N=0 

YLMsYLIM) 

CMxCOSI ACR4YLM) 

80  N=N41 

IF  (N-N9)  81.81*75 

81  L=L4l 
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K=0 
N0O=0 
XLN=XLI N) 

IF  (KQ-1)  83.S2.83 

82  NOOsO 

OO  182  |K=1,S 

182  OVR(IK)-0.0 
GO  TO  84 

83  OO  183  IK=6.18 

183  OVRllK)sO.O 
«=■! 

ONH=M 
NOD=  1 
K=0 

OHM=Hl.<Ll 
OUH=UL(L) •« 

0VH  = V1_<L>4W 
OTH  = n.<L)  •• 

B4  K=K«-I 

IF  (K-KS)  85.85.90 

384  I=14KS-K 
GO  TO  90 

85  I=I+l 

IF  (NOO-KQ5)  86.384.384 

86  J=IS< I ) 

IF  (J)  84.  84.87 

87  IF  lABSIUSUn.GT. 100.0)  GO  TO  84 
XSJ=XS<  J» 

YSJ=VS( J) 

LlSJ  = US«  J) 

VSJ=VS<  J) 

HSJ=HS( J) 

TSJ=TSC  J) 

DYS=YSJ-rLM 

DXS=  (XLN-XSJ)4CM 

0YS2=0YS4DYS40XS40XS 

0X52=0. 540YS2 

IF  IIOS.EQ.O)  GO  TO  385 

USK=USJ 

VSK=VSJ 

IF  IKQ.EQ.l)  GO  TO  388 
USKsUSJAULI L ) 

VSK=VSJ4VL(L) 

388  OXS1SUSK4USK4VSK4VSK40.01 
0XS2=( USK4DYS-VSK4  0XS  ) 
0X52=0X52*0X52/0X51 

385  *=  C2/t 0Y5240X52*ALPH4  C2) 
NOO=NOO«-1 

0NH=0NM4I« 

T5J=T5 J*H 
H5 J=H5 J*W 
U5J=USJ*M 
V5 J=V5 J*W 
0HH=0HM4H5 J 
0UF=0UM4U5J 
0VH=0VH4VS3 
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c 

c 

c 

c 

c 


OTK=DTH+TSJ 
IF  (KU-I)  84,84«a9 

89  DYHsDYH+DVS** 

OXHsOXH«-OXS*N 
OXYH=D  XYH  + OX  S*DYS*  4 
0XXH=0XXH40XS*DXS*» 

OY  yH=nYYH-*’OYS*OYS*W 
OXHH=DXHH+HS J4DXS 
DYHH=DYHH+MSJ*DYS 
OXTH=D  XTM+TS J4DXS 
OYTH=DVTH-'.  VSJ40YS 
OXUHsOXUH^US J*OXS 
DYUH=OYUH+US J4DVS 
OXVH=DXVH+VS J*OXS 
OYVH=DY  VH-fVS  J*OYS 
GO  TO  84 

90  CONTINUE 

IF  (KQ-1)  94.91*94 

91  IF  lNOD-21  80*92*92 

92  IF  (ONHI  80.80*93 

93  0NH=1.0/0NH 
HL(L)=0HH4DNH 
Ul.lL)=OUH*ON« 

VL IL I-DVH4DNH 
TL<L)=DTH*ONH 
GO  TO  80 

94  IF  (NOD-3)  80*80*95 

95  D=0YH*DYH-DNH4DYYH 
E=DXYH4DVM-0XH*DYVH 
4=0XH*0YH-DNH*0XYH 
B=OXXM*OYH-OXH*OXVH 
BOAE=S4D-A*E 

IF  (BOAE)  97*80*97 

97  BI=1.0/BOAE 

C =DXHH*OYH-OMH*DXYH 
F =DYHH*DYH-DHH*DYVH 
CT=OXTH4DYH-OTH*OXYH 
FT=DYTH*OVM-OTH*DYYM 
CU»0XUH40YH-0UH*0XYH 
FU=OYUH*DYH- 0UH40YYH 
CV=0XVH*DYH-DVH4DXYH 
FV=0YVM4DYH-OVH4OVYH 
HL (L)=(B4F-C*E)4BI 
UL(L)= ( B4FU-CU*E>4BI 
VL (L  >*<B4FV-CV4E»*BI 
TL(L)  = < B4FT-CT*E>*8I 
GO  TO  80 
100  CONTINUE 
RETURN 

THIS  SUBSECTION  DETERMINES  A SET  OF  KS  CLOSEST  OBSERVATIONS  FOR 

EACH  GRID  POINT.  THE  OBSERVATIONS  INOENTIFVING  INDICIES  ARE  STORED 

IN  ARRAY  IS  IIS  DIMENSIONS  MUST  BE  GREATER  OR  LESS  THAN  LRH.C4KS) 


18  ACR=3. 1416/180.0 
OLCK=RMAX4RMAX 
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JCTa  50 
JS<  JCT)=-I 
DSC  JCT  ) = OLCK 
JC2=JCT-2 
M=0 
L=0 
1 = 0 

20  M=M+l 
N=0 

YLM=VL  CM  > 

CM=COSC  ACR*tC.M) 

25  N=N«’l 
L=L>» 

J=0 
JT*J  J 
XLN=XL  CN) 

XLCK=OLCK 

29  J=0 
JC  = 0 

30  J=J+l 
DyS=YSC  J »-YLM 
DXS=CXLN-XSC Jl )*CM 
DLS=OYS*OVS+DXS*DXS 

IF  COLS-XLCK)  32,33«33 

32  JC=JC+i 

IF  CJC.LT.JCTI  GO  TO  31 

XLCK=XLCK-I0.0 

GO  TO  29 

31  OSCJC)=DL.S 
JSC  JC)=J 
IYS=0 
IXS=l 

IF  (DXS.LT.O)  IXS=2 
IF  COYS.LT.O)  IYS=2 
IQC JC)=IYS+IXS 

33  IF  <J-JTI  30.34*34 

34  KP=0 
K=0 
JQ=1 

I SM=1 
JX  = JC 

IF  (KSD.LT.l)  ISW=2 

35  K=K4-1 

S IL=DLCK 
J=0 
I = I41 
KP-KP-f  1 
JC= JCT 

40  J=J*l 

IF  CJ.LE.JX)  GO  TO  C4l.42}.lS«f 
GO  TO  45 

41  IF  1 IQC J I.NE.JQ)  GO  TO  40 

42  IF  COSfJI-SIL)  44.40.40 
44  SII.  = OSCJ) 

JC=  J 

GO  TO  40 
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45 


49 

50 
55 
60 


OS( JC) sOLCK 
1S( I ) = JSC  JC) 

J0=JQ41 

IF  CJQ.GT.4)  JQs| 

IF  (KP.EQ«KSW)  |SII  = 2 
IF  CJC.LT.JCT}  GO  TO  49 
IF  (KP.GT.KS«)  GO  TO  49 
K=K-I 
I = I-l 
GO  TO  35 

IF  (K-KSI  35*50*50 

IF  (N-N9)  25*55.55 

IF  (M-M91  20*60*60 

CONTINUE 

HE TURN 

END 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUbROUTINF.  CHECK  ( J J . YS  t XS  . VS  t US  , VL  . XL  t V N , UN  ) 

THI SSUbROUTINE  CHECKS  THE  CSSERVEO  WIND  VALUES  HY  COMPARING  THEM 
TO  ANALYZED  A I NO  VALUES  I THE  U AND  V OF  THE  INCONSISTENT  DATA  ARE 
SETAT  «i99.y  ). 

JJ  = NUMUER  OF  WIND  DATA 

YS.XS  = LATITUDE  AND  LONGITUDE  OF  WIND  DATA  lOFGI 

US«VS  = U AND  V COMPONENTS  OF  WIND  DATA 

YL.XL  = LATITUDE  AND  LONGITUDE  OF  ROWS  AND  COLUMNS 

UN.VN  = GRID  POINT  U AND  V WIND  COMPONENTS 

CRT=  CRITICAL  VALUE  USED  IN  TESTING  WIND  DATA 

WIM=  WEIGHT  GIVEN  TO  A WIND  IN  ANALYZING  A VALUE  AT  ITS  LOCATION 

AND  WEIGHT  GIVEN  TO  NEAREST  GRID  POINT  VALUE  I NOW  CALLED  SIM  AND 
GIM) 

LOGICAL  DEBUG 

DIMENSION  JSTI  26).0ST<  26)«WTS(  25)tOVR(20) 

DIMENSION  XSI  D.YSI  II.VSI  DtUSI  1) 

DIMENSION  XL!  1 ) « YL ( I).UN(  I ) • VN ( 1) 

COMMON/CIS/  ISIIOOO) 

COMMON/CCHK/  CRT. SI M, G I M* I WND 

COMMON /CRTS/  KS.Wl»C2.RMAX»KSS5»IDS.KSW,ALPH 
CUMMON/CGO/  M9.N9.  I 9 . M8 . N8 . YB . XB. DD 
EQUIVALENCE  ( USK . US J I . I VSK . VS J ) 

equivalence  (DVRI I ) ,ONH) . <DVRI2 ) .DHH) , ( DVR ( 3 ) . DUH ) . (OVR(A ) . DVH ) . 

2 (DVRI 5 ).DTH) . ( DVR (6) .OXH> . (DVRI 7 ) .DYH) . ( DVR( 8» .DXYH) . 

3 (DVR(9 ) ,OXXH) , (DVRI  I 0) .DYYH) . IDVR( 1 1 ) .OXHH) . ( DVRI  12) . DXUH) , 

4 (DVRI  I 3) .DXVH) . (DVRI  14) .DYHH) . (DVRI I 5 ) . DYUH ) . I DVRI 1 6 ) ,OYVH) , 
b (DVRI 17).OXTH). (DVRI 18) »OYTH) 

DATA  DEBUG/. FALSE./ 


C 

C 

c 


c 

c 

c 

c 


4 FQRMAT(2l X. . 13H  DELETED  DATA/) 

5 FORMAT  I49H  LAT(OEG)  LQNIDEG)  UIKTS)  VIKTS)  TEST./) 

6 F0RMATUX.I0F10.2) 

PRINT  4 

PRINT  5 

BASIC  COMPUTATIONS 
JT  = JJ 

ACR=3. I4I6/I80.0 
YOI=l .O/DD 
XDI=1.0/OO 
KQ5=KSS542 

IF  (IWND.GT.O)  KQ5=KQS41 

SPHI=ALPH4ALPH 

CRS=1 .254CRT 

ANALYSIS  OF  WINDS  AT  MEASUREMENT  LOCATIONS  (SIMILAR  TO  METHOD  USED 
IN  SUBROUTINE  MESH). 


20 


KG=0 
JX=0 
LX  = 0 
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CRW-CR5 

IF  (KG.EQ.2)  CRW=CRT 
80  LXsLX^l 

IF  (LX-JT)  82.82*100 
82  L=1.X 

IF  (KG.NE.O)  L-  JST(LX) 

IF  ( ABSIUSIL) 1 .GT. 100.0)  GO  TO  80 
YLM=VS(L) 

XLN=XSIL) 

M=<  YLI  1 )-rLM)*yoi-n  .5 
N=(XLN-XL<1  ))*XOI  + i.5 
IF  CM.LT.l)  M=1 
IF  (M.GT.H9)  M=M9 
IF  IN.LT.I)  N=1 
IF  (N.GT.N9)  N=N9 
l.I=N9*(  H-1  ) ♦N 
I=ILI-1 )*KS 
CM=COS< ACRFYLM) 

NOO^O 
K = 0 

00  182  1K=1«18 

182  OVR(IK)=0.0 

IF  IIWNO.LE.O)  GO  TO  84 
XSJ=XL<  N) 

YSJ=YCCM) 

US  J-=UN(LI  ) 

VSJ  = VN(LI  ) 

DYS=YSJ-YLM 
OXS=( XSJ-XLN ) ACM 
W=GIM 
GO  TO  89 
B4  K=K41 

IF  (K-KS)  85.85.90 

85  I=I41 

IF  (NOO-KQ5)  86.84.84 

86  J= IS(  I ) 

IF  <J)  84.84.87 

87  IF  f ABSIUSI  JD.GT.IOO.O)  GO  TO  84 
USJsUSlJ  ) 

VSJ  = VS<  J) 

XSJsXSI  J ) 

YSJ  = YS<  J ) 

DYS=YSJ-YLM 
DXS=<  XSJ-XLN)4CM 
W=SIM 

IF  (J.EQ.L)  GO  TO  389 
DXS2=0.0 

DYS2=DYS*OVS4DXS*DXS 
IF  IIOS.LT.l)  GO  TO  88 
OXS1SUSK4USK4VSK4VSK40.01 
DXS2=( USK4DYS-VSK4DXS) 
OXS2=OXS2«OXS2/OXSl 

88  W=C2/(  0YS24-0XS24SPH14C2) 

389  IF  (KG.NE.O)  WsW*MTS(J) 

89  NaD=N0041 
DNHSONH4W 
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US J=USJ*W 

vs J=VSJ*W 

OUH=DUH«-USJ 

OVK=OVH+VSJ 

OVH=DYH>OVS*X 

OXH=OXH*DXS*W 

OXVK=OXYH*OXS*DYS*W 

OXXH=OXXH+OXS*OXS*l« 

L)YYK=DYYH-*-DYS»OYS*ll( 

OXUH-=DXUH  + USJ*DXS 
DYUH=OYUH+USJ*DYS 
DXVH=OXVH+VSJ*OXS 
OYVH=DYVH*VS J*OYS 
GU  TO  84 
QO  CONTINUE 

Q4  IF  (NUO-3)  80.80.95 

95  0=DYH*0YH-0NH*DYYH 
E=OXYH*DYH-DXH40YYH 
A=DXH*DYH-ONH*DXYH 
a=OXXH*OYM-DXH*OXyH 
bOAE=B*I>-A4E 
IF  (BOAE)  97.80.97 
97  BI=l.0/8OAE 

XU=0XUH40YH-DUH4DXYH 
YU=OYUH*OYH-OUH*OYYH 
XV=DXVH*OYH-OVH*DXYH 
YV=DYVH*OYH-OVH*DYYH 
ULL  =C B*YU-XU*EI4BI 
VLL  =( B»YV-XV*E ) ♦BI 
C 

C CHECK  FOR  FOR  INCONSISTENCY  BETWEEN  ANALYZED  AND  MEASURED  WINDS 
C 

ALTHl=  ULL*ULL+VLL*VLL 

ALTH2=  US<L)*USIL)*VS(L)*VS(L) 

aLTH=ALTHl 

IF(ALTH2  .GT.  ALTHl)  aLTH=  ALTH2 
PERPl=( ULL*US(L ) ♦VLLWVSIL) I/BLTH 
IF  (DEBUG) 

2PRINT  6.  YS(L)  .XS(L)  .USID.VSID.PERPl  .ULL.  VLL.CRW.CR5.CRT 
DLS=PERP1 
WTS(L)=1 .0 

IF  (OLS.GT.CR5)  GO  TO  80 
WTS(L) =0.5 

IF  (DLS.GT.CRW)  GO  TO  80 
IF  (KG.EQ.2)  GO  TO  99 
WTS(L)=SI H 

IF  (JX.GT.25)  GO  TO  100 

JX  = JX4'1 

JST(JX)  =L 

DST(JX)  = PERPl 

GO  TO  80 

99  PRINT  6.  YS(L)  .XS(L)  .USID.VSIL). PERPl 
US(L)  =-999.9 
VS<L)  =-999.9 
GO  TO  80 
too  CONTINUE 
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C ORDER  SUSPECT  DATA  ACCORDING  TO  HIGHEST  PERPl  VALUES  AND  REPEAT  AHOVE 
C 

IF  (JX.LE.O)  GO  TO  300 
JT=JX 

IF  (KG.EQ.O)  GO  TO  202 
IF  IKG.NE.l)  GO  TO  300 
DO  200  J1=1*JT 
J3  = JI 

DO  150  J2=Jlt3T 

IF  (0STIJ2)  .LT.  0ST(J3n  J3=J2 
150  CONTINUE 

IST=JST( J1  ) 

XST=OST{ Jl  ) 

JSTIJl »=JST( J3» 

OSTI Jl  ) = OST< J3  ) 

JSTI J3)=I ST 
OSTI J3  > = XST 
200  CONTINUE 
202  KG=KGF1 
GO  TO  20 
300  CONTINUE 
Rf  TURN 
END 


SUBROUT INE  K10(l V* I 0*  I B* VOR • O 1 V . BAL • U « V * C . T ) 


C 

C THIS  SUBROUTINE  COMPUTES  FIELDS  OF  VOR TI C I TY , DI VERGENCE . AND  THE 
C HALANCe  HEIGHT  TERM. 

. C 

C IV. ID. IB  = COMPUTES  VORT I C IT Y. D I VERGE NCE  AND  BALANCE  TERM  WHEN 
C IV.ID  AND  IB  ARE  SET  GT  0 

C VOR.DIV.BAL  = VORT I C I TY , D I VERGENCE  AND  BALANCE  TERM  VALUES 
C U.V  = U AND  V WIND  COMPONENTS 

C C.T  = LATITUDINAL  COSINE  AND  TANGENT  VALUES 
C 

DIMENSION  VORin.OIVil  I.BAL  (D.UIII.VIII.CI  1).T(  1) 
COMMON/CGD/  M9.N9.I9.M6.N8.TB.XB.D0 

DATA  ACR.G.R.TDEL  /O . 0 17A533 . 9 .8062. 6 . 37  IE 6 , I . ASBAE- ♦/ 

ERI=l.O/R 

GI= I .O/G 

OY=DO*R*ACR 

E DYI=1.0/f DYADY) 

I  VB=0 

IF  ( I V.GT. O.OR.IB.GT.O ) I VB= 1 
t DO  50  M-2.M8 

CM=C(M) 

DYI S=DYI/CM 

TNGR=T (M)«ER1 

IF  (IB.LT.I)  GO  TO  10 

CM2=0.5*GI ACM 

HC2=DV ACMADYACM 

BFM  =GI ATDELAERIACM 

CFM  =GIATDELACMAT(M) 

10  DO  40  N=2.Ne 
I=<M-1 IAN9+N 
IF  1 ID.LT. II  GOTO  20 
DUE=OYISA<Ul !♦ l)-U( 1-1)1 
OVN^OYlAlVf I-N9)-V<  IAN9)  ) 

01 V( I l=DUE+DVN-V( I lATNGR 
20  IF  (IVB.LT.l)  GO  TO  AO 
OVE=DYISA(VI I+l )-V( I-l I) 

DUN=DY1 A( U( I-N9)-U( I AN9) I 
VORI  = DVE-DUNAUI 1 lATNGR 
IF  (IV.GT.O)  V0R(1)=V0RI 
IF  ( IB.GT  .0  I 

IBALII)=  tVORIACFM  -UIIIABFM  )AHC2 

2 -CM2A(  (V1IA1)-V<1-1)|A(U(  I-N9)-U<  I AN9 I I 

3 -(U(  1*1  )-U<  I-l  I )A(  V(  I-N9  »-V(  IAN9)  I ) 

40  CONTINUE 

SO  CONTINUE 
RETURN 
END 
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SUBROUTINE  ALTERSt LGM «CT V «F AC • TNG *05 • UN • VN« VV .OV ) 


C 

C THIS  SUBROUTINE  CHANGES  THE  INITIAL  WIND  FIELD  INTO  A NONOIVERGENT 
C FIELD  THAT  FITS  THE  SPECIFIED  VORTICITY  AND  DIVERGENCE  FIELDS. 

C 

C SPHERICAL  GRID  VERSION. 

C 

c LGM  = MAXIMUM  NUMBER  OF  ITERATIONS. 

C CTV  = ERROR  TOLERANCE  IN  VORTICITY  AND  DIVERGENCE  FIELDS  (SEC-1). 

C FAC  = MIND  ADJUSTMENT  FACTOR. 

C CS.TNG  = ARRAYS  CONTAINING  LATITUDINAL  COSINES  AND  TANGENTS. 

C UN.VN  = ACTUAL  WIND  COMPONENTS  - INITIAL  VALUES  (M  SEC-t). 

C UN.VN  3 NONOIVERGENT  MIND  COMPONENTS  - FINAL  VALUES  (M  SEC- 1 ) . 

C VV  = VORTICITY  (SEC-1). 

C OV  = DIVERGENCE  (SEC-1  ) 

C 

LOGICAL  DEBUG 

DIMENSION  TNG( 1 ) .CS ( 1 I .UN(  1 }.VN( 1 ). VV( 1 ).OV(  1 I 
COMMON/CGD/  M9 • N9. I 9. M8 . NB. YB. XB.DO 
EQUIVALENCE  (LR .M9  ) • ( LC. N9  1 
DATA  DEBUG/. false./ 

C 

C PRELIMENARY  COMPUTATIONS. 

C 

R=6.371E6 

0Y=00*R*3. 1416/180.0 
LG=0 

0Y2=0Y+DY 
0YRA5=DV*FAC 
ER Isl .0/R 
OYI=l .0/DY2 
LRl=LR-l 
LCl=LC-l 

XRCI=l. 0/(LH*CC) 

COMPUTATION  OF  AVERAGE  U AND  V VALUES. 

SUMlsO.O 
SUM2=0.0 
1=  1 

DO  8 JsI.LR 
DO  8 Ksl.LC 
SUMls  SUM  I ♦ UNI  I) 

SUM2=  SUM2  AVNI  I I 
8 I=l*l 

SUM1=SUM1*XRCI 
SUM2sSUM2*XRCI 

BEGINNING  OF  ITERATION  LOOP. 

10  LG  =LG41 
OIFDa^O.O 
01FV=0.0 

BEGINNING  OF  GRID  POINT  LOOP. 


j 
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DO  500  J=2,URl 
TNG637=TNGC  J » *ER I 

DrisEc=oyi/cst  j) 

DYRACS=0YHA5*CSt  Jl 
DO  500  K=2«LCl 
I=( J-1 ) AtC+K 
IPi=I+l 

IP9=I-*-N9 

IM9=I-N9 

C 

C ALTER  WINDS  TO  FIT  A ZERO  DIVERGENCE. 

C 

DUE  = DYISEC*lUN(  IPD-UNI  IMD) 

DVN  = DYI*(VN(  IM9)-VNI  IP9  ) ) 

DIJK  =DUE+OVN-VN(  I)*TNG637 
DIF=DV(  I )-DIJK 

IF  ( ABSIDIFI.GT.DIFO)  DIFD= ABSI DI F ) 
CUJK  =OYRACS*OIF 
CVJK  =0YRA5*0IF 


VNI 

IP9> 

= VNI 

I P9) 

- 

CVJK 

VNI 

IM9) 

= VNI 

IM9) 

4- 

CVJK 

UNI 

IPI  ) 

= UNI 

IPI  I 

♦ 

CUJK 

UNI 

IMl  ) 

= UNI 

IMl  I 

- 

CUJK 

C ALTER  WINDS  TO  FIT  THE  SPECIFIED  VORTICITY  VALUE. 
C 

DVE  = DY  ISEC*I  VNI  IPI)--VN*  IMl)) 

DUN  = DYI*(UN(  IM9)-UNI  IP9>) 

VOJK  =DVE-DUN+UNI  I)*TNG637 
VIF=VVI  II-VOJK 

IF  ( ABS( V IF  I.GT.D IFV)  O I F V= ABSI VI F ) 

CVJK  =OYRACS*VIF 
CUJK  =DYRAS*VIF 


UNI 

IP9) 

= UNI 

IP9) 

♦ 

CUJK 

UNI 

I M9} 

= UNI 

IM9) 

- 

CUJK 

VNI 

IPI  1 

= VNI 

IPII 

* 

CVJK 

VNI 

1 Ml  ) 

•=  VNI 

I Ml  I 

- 

CVJK 

500  CONTINUE 
C 

C CHECK  OF  ITERATION  NUMBER  AND  WAX I MUM  DIFFERENCE  BETWEEN  THE  SPECIFIED 
C VORTICITIES  AND  DIVERGENCES  AND  THOSE  COMPUTED  FROM  THE  WIND  VALUES. 
C 

IF  (LG.GE.LGM)  GO  TO  600 

IF  (DIFO.GT  .CT V.OR .DIFV.GT .CTV)  GO  TO  10 
bOO  CONTINUE 
C 

C ADJUSTMENT  OF  U AND  V FIELDS  TO  AVERAGE  VALUES  OF  INITIAL  FIELDS. 

C 

SUM3=0.0 

SUM4=0.0 

I=I 

DO  810  J=I«LR 
DO  eio  K-sI.lC 
SUM3=  SUM3  «■  UNI  II 
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SUM4=  SUM4  ♦VNC  I) 

810  I=I*1 

SUM3=SUM34XHCI 
SUM4=SUM4«XRCI 
SUM5=SUMI-SUM3 
I SUM6=SUM2-SUM4 

I=l 

00  820  Jsl.LR 
OO  820  K=1.LC 
UNI  1 )=  UNI  I )4SUM5 
VN(  I)=  VNI  Il-t-SUM6 

820  I=I+l 
C 

PRINTOUT  OF  MAXIMUM  DIFFERENCES  AND  COMPUTER  TIME. 

IF  I DEBUG) 

2PRINT  2000.LG.DIFV.DIFD.RT2 

2000  FORMAT ( 2X, "NO  OF  ITERATIONS  = ".I3.5X,"MAX  VORT  ERROR  E 1 1 .2 . 5X , 

2"MAX  OlV  ERROR  =" . El  I . 2. 5X . "CP  TIME  =".E1I.2) 

C 

C END  OF  SUBROUTINE  ALTERS. 

C 

RETURN 
END 
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SUBROUTINE  BALHGT (NR. NC»  OL AT . I BT . N IT , «C . EPS . AV E. CF . TF . HU. XOR) 

THISSUBROUT INE  SOLVES  THE  POISSON  EQ.  (LAPLACIAN  OF  HU=XOR  I FOR  THE 
SCALAR  FIELD  ANAL  USING  THE  ALTERNATING  DIRECTION  IMPLICIT  (ADM 
METHOD  OF  RELAXATION  ( YOUNG . D. M. . SUR VEY  OF  NUMERICAL  ANALYSIS* 
MCGRAW-HILL  .1962 . EDITED  BY  J.TODD).  A FIRST  GUESS  FIELD  SHOULD  BE 
SUPPLIED  FOR  HL.  HOWEVER. AN  EQUIVALENT  RESULT  CAN  BE  OBTAINED  WITH 
ADDITIONAL  ITERATIONS. 

NR.NC  = TOTAL  NUMBER  OF  ROWS  AND  COLUMNS  IN  GRID. 

OLAT  = UNIFORM  SPHERICAL  MESH  SIZE  IN  DEGREES. 

IBT  = INDICATOR  SPECIFYING  TYPE  OF  BOUNDARY  CONDITIONS. 

(IBT.LE.OI  MEANS  FIXED  BOUNDARY  CONDITIONS. 

(IBT.GT.O)  WEANS  NORMAL  DERIVATIVE  BOUNDARY  CONDITIONS. 

NIT  = NUMBER  OF  ITERATIONS  TO  BE  PERFORMED. 

WC  = RELAXATION  FACTOR. 

EPS  = RELAXATION  TOLERANCE  -IF  THE  CHANGES  IN  HL  DURING  THE  ROW 
IMPROVEMENTS  OF  THE  NTF  ITERATION  ARE  LESS  THEN  EPS. 

THE  RELAXATION  IS  STOPPED  AFTER  THE  (N+IITH  ITERATION. 

AVE  = AVERAGE  VALUE  SPECIFIED  FOR  ML  - THE  FIELD  IS  ADJUSTED 
TO  THIS  AVERAGE  VALUE  IF  (IBT.LE.OI 

CF.TF  = TWO  ARRAYS  CONTAINING  LATITUDINAL  SINE  AND  TANGENT  VALUES. 
HL  = THE  SCALAR  FIELD  TO  BE  DETERMINED  IM2SEC-I  OR  M) . 

XOR  = THE  FORCING  FUNCTION  (M2SEC-I  OR  M). 

LOGICAL  DEBUG 

DIMENSION  TY(  10) .BB(  10) *B1 ( 1 0 ) *B2 ( I 0 ) . E I 1 0 ) *F(  10 ) 

DIMENSION  HL ( 1 ). XOR( 1 >.CF(  I )*TF(  1 ) 

DIMENSION  80W(IO)*BOE(IO) .BDNI 10) . BOS ( 10) 

PRELIMENARY  COMPUTATIONS. 

DATA  MI. Nl.ACR/2. 2*0. 0174533/ 

DATA  DEBUG/. FALSE./ 

M9=NR-1 
N9=NC-I 
I 9SNR4NC 
I9M= I9-N9 
ICS=I 
IRS=NC 

IF  (IBT.LE.OI  GO  TO  100 
DO  62  MSM1.M9 
I = ( M- 1 I *IRS^1 
BOW( MI=HL ( I )-ML( I ♦ICS ) 

I=I4.N9*ICS 

62  BDE(M)sHL(  I )-HL(I-ICSI 
DO  64  N=Nl .N9 
I=(N-I >*ICS4 1 
BON(N)=HL ( I )-HL(  I ♦IRS  > 

I=I4M9*I RS 

64  B0S(N)  = HL(  I l-HL<I-IRS) 

G=9. 80  62 
TOEL= I .4S84E-4 
R=6.371E6 
CR=3. 1416/180. 0 
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CON=l .0 
CIR^O.O 
V0R=0.0 
00  70 

CSI=i.0/(CF4M) •CF ( M ) • TF I M ) ) 
C1R=CIR'»'BD«(M)*CSI 
70  CIH-=ClR<-aOE|M»*CSI 

51=2.0/ (TF (Ml )»TF(M1-1 ») 
59=2. 0/( TF(M9»»TF( M9^1 )} 

DO  75  N=N1.N9 
C1R=CIR^B0N(N)«51 
75  CIR=CIR^BD5CN» *59 
C1R=CIR*C0N 
00  80  M=M1.M9 
FCM=1.0/CF(M) 

00  80  N=N1.N9 

1 = ( M-1  )*NC'fN 

80  VOR=VOR«'XOR(11*FCM 

DIF=( VOR-CIHI *0.5/1 M9+M9-2I 
DIF=OIF/CON 
OC  85  M=M1.M9 
C5I=CF(M) *CF(M)*TF<M> 

BOM( M)  = BD«f  M1*01F*CS1 
85  BOE(M)=BOE(M}4-OIF*CSl 

51={ TF( Ml ) ♦TF( Ml-1 } ) *0.5 
59=(TF<M9 )*TF( M9*l ) )*0.5 
00  90  N=N1.N9 
BDN(N)  = BDN  (N  )'*'0(F*51 

90  B05(N)=B0S1  M-*'0IF*59 
CIR=0.0 

OC  91  M=M1.M9 

C5I=l.0/(CF(M)*CF(M)*TF<M»l 

CIR=CIR*BOW<Ml*CSI 

91  CIR=C1R4^BDE(M)  *C5I 

51  = 2.0/(TF(M1  )+TF(Ml-l  )) 
59=2.0/4 TF«M9)*TF I M9^ I )) 

OO  92  N=N1.N9 
C(R=CIR^B0N4N>*S1 

92  CIR=CXR*BD5(N)*S9 
C1R=C1R*C0N 

100  0H=0LAT«ACR*0«S 
00  150  M=M1«M9 
BB(M)3CF4M)*CF4M) 

150  TV4M)=0H*CF(M)*CF4M>*TF<M| 

NI  = 0 
MP=NIT 

180  «CI=1.0/HC 
I)C1  = 1 .0-MC 
OO  190  M=M1.M9 
B1  (M)=(BB(  M}4-BB(  M)*2.0  )*WCI 
190  B24M)=B1 (M)*«C1 
C 

C BEGINNING  OF  ADI  ITERATION  LOOP. 
C 

200  NI=NI>1 
C 
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successive  how  improvement. 

RLI  = 0.  0 

OO  21 S M=MI«M<) 

1=<M-1 )*IRS«I 
BBM^BBIM) 

BIM-Bl IM) 

B2M=B2( M) 

TYM=TYf M) 

BPT=BBM+TYM 
aMT=BBM-TYM 

IF  IlBT.LT.l)  GO  TO  204 
El  1 )=  1.0 

FI  n=  BDWIM) 

GO  TO  206 
204  Ell)-  0.0 

FI  1 )=  HLI  1 ) 

206  1=H-ICS 

00  210  N=N1«N9 
208  e(N)=  1.0/ieiM-EIN-l )) 

F|N>  = lHL(  I+1RS)*8PT»HCII-IRS)4BMT*B2M*HH  I ) -XOH I I ) +F ( N- 1 ) ) WEI N ) 
210  I=I+ICS 
NB=N9+1 

IF  (IBT.GT.O)  HLI I )=l BDEI M) ♦FIN9) )/( 1 .0-ei N9) ) 

00  213  N=N1.N9 
I=I-ICS 
N&=NB-  1 

OLI=HLll )-E<NB)4HL( I*1CS)-FINB) 

ALI=ABS( OLl ) 

IF  lALl.GT.HLI)  RL1=AL1 

213  HLI 1 )=HLI I )-OLl 

IF  IIBT.GT.O)  HLI 1-ICS)=E|NB- 1I4MLI I )4FlNB-l I 

215  CONTINUE 
C 

C SUCCESSIVE  COLUMN  IMPROVEMENT. 

C 

00  240  N=N1.N9 
l=|N-l )*1CS41 
IF  IlBT.LT.l)  GO  TO  214 
Ell)=  1.0 
FIDS  BONIN) 

GO  TO  216 

214  El  1 )s  0.0 
Fll)=  HLI I I 

216  Isi  + IRS 

00  230  MSM1.M9 
218  TYMsTYIM) 

BMTsBBI M)-TVM 
EIsl.O/IBl I M)-BMT*EIM-1 ) ) 
eiM)=  IBBIM)*TYM)*EI 
FI M)sFT 

FT  s IHLI  I4-1CSD-HL<  I-tCS)4B2  IM)  4HLI  1 )-XOR|  1 )4^F  IM-l  ) 4BMT  ) 4E  I 

230  Is  DIRS 
MBsMOl- 1 

IF  IIBT.GT.O)  HLI I)s|B0SIN)^FIM9))/I 1.0-EIM9)) 

00  235  MsMI  .mo 
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I=I-IRS 
MU=MB-1 
23S  HL(  I ) = 

IF  (IBT.OT.O)  Hl.(  I-IRSME(  MB-1  I^HL  (1  I 

240  CONTINUE 
C 

C CHECK  OF  THE  TOLERANCE 

C ANO  THE  ITERATION  NUMBER  I THE  RETURN  IS  NORMALLY  MADE  TO  200*  BUT  AT 

C THE  NEXT  TO  LAST  ITERATION  WC  IS  SET  TO  0 ANO  RETURN  IS  MADE  TO  160) 

C 

IF  (NI41 .GT.NIP)  GO  TO  285 

IF  I NIAl .LT.NIP.ANO.RL  I.GT.EPS » GO  TO  180 
MC3:1.0 
NIP*NI+1 
GO  TO  180 
285  CONTINUE 
C 

C COMPUTATION  OF  RESIDUALS. 

C 

SUM=0.0 
RRS=0.  0 

saR=o.o 

OO  275  M=M1«M9 
1=(M-I)4IRS42 
BBM=BBf M) 

BB2=BBM'l-BBM4-2.0 
TVM=TY<M ) 

FCM=1. 0/CFI N) 

8PT=BBM+TVM 
aMT=BBM-TYM 
OO  270  N=NI.N9 

RRSS=CHL(  f 4ICSIFHL  { I-ICS  )-»HL(  I -HRS  i *SPT*HL(  l~tRS}  *BI4T 
2-8B2*HL( 1 )-XOR( I ) } •FCM 
RRS=RRSARRSS«RRSS 
SUM-SUMAl .0 

SOR=SOR4ABS( XORI I ) 4FCM) 

270  I=I4ICS 
275  CONTINUE 

RR S=SQRT(RRS/SUM) 

SOR-SOR/SUM 
PCT=1 00. 04RRS/SOR 
C 

C PRINTOUT  OF  RESIDUAL  ANO  COMPUTER  TIME  INFORMATION. 

C 

IF  fOEBUG) 

2PRINT  500.SOR.RRS.PCT.  N1 
500  FORMAT (69H  AVE  ABS(FOR)  RESIDUAL  PERCENT 

I NO  OF  ITER./ IX* E12.4.3E14.A. 18 ) 

C 

C ADJUSTMENT  OF  SCALAR  FIELD  HL  (ANAL)  TO  THE  SPECIFIED  AVERAGE  VALUE. 

C 

IF  I 18T.lt.  I)  GO  TO  400 

HLI=0.0 

SUM=:0.  0 

OO  390  MSM1.M9 
I=C  M-l  )*1RS41 
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00  390  N=N1*N9 
I=I4^ICS 
SUM^SUM-t-l  .0 
390  HH=HLI+HL<1) 

HLI=MCl/S«JM 
HLI  = AVe-Hl_I 
00  395  M=1.NR 
I=<M-i  )*IRS«-1 
OO  395  N=».NC 
HI-(  1 ) = HL<  I )»HLI 
395  I=I^ICS 

HL( I )=(HL<2)+HL« 1 ♦NC) 1*0.5 
HL<NC)  = C ML( NC- 1 )*HL* NC+NC) » *0.5 
HL ( I9MI=( HL( 19M+I ) ♦HL ( I9M-NC) )*0.5 
HL(  I9)  = (Hl.l  19-1  14-M4.I  I 9-NC  ) ) *0.5 
♦00  CONTINUE 
C 

RETURN 

ENO 
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SUBHOUT INE  AO VEC( DT . US « VS *0  1 .02 •CM) 

C 

C THIS  SUBROUTINE  AOVECTS  THE  FIELD  01  WITH  A SPECIFIED  WIND. 
C 

C DT  = TIME  STEP  OVER  WHICH  AOVECTION  IS  MADE 
C US* VS  = COMPONENTS  OF  WIND  USED  TO  AOVECT  01  FIELD 
C 01  = FIELD  THAT  IS  ADVECTED 
C 02  = DUMMY  FIELD 
C CM  = COSINES  OF  ROW  LATITUDES 
C 

DIMENSION  OK  1)«02(  1)«CM(  1 ) 

COMMON/CGO/  M9.N9. I 9. MB.NB. VB* XB . DO 

XlD-1 . 0 

YI0=XID*N9 

TT=OT 

UT  = US 

VT  = VS 

IT=i 

DO  20  IR=  1*  10 

IF  (UT.LT.l .O.AND.VT.LT.l  .0  ) GO  TO  22 
UT=0.5WUT 
VT=0.5*VT 
TT=0.5WTT 
20  IT=IT+1 
22  CON  =ABS{TT) 

DO  200  1R  = 1 • I T 
1=  1 

DO  100  M=1.M9 
CMI=1. O/CMIM) 

DO  100  N)=1«N9 
011=0111) 

UA=ABS(UT)*CMI 

VA=A8S(VT) 

IF  (DT.LT.O.O)  GO  TO  80 
1U=N-S1GN{XI0«  UT) 

1 V=l+SIGNCy 10. VT) 

GO  TO  90 

no  IU=N'I-S1GN(X10.UT) 

I V=I-S1GN<  VIO. VT) 

90  01V=0.0 

01 U=0.0 

IF  ( IV.GT.O.ANO.IV.LE.  19)  01V=Q1(IV) 

JU=IU*<M-1 )*N9 

IF  ( 1U.GT.O.ANO.IU.LE.N9)  aiU=QllJU) 

02(  I ) = 01  I-(  (01  I-OlU  )*UA4-(Ql(-0  IV  )*VA)«>CON 

100  1=14-1 

DO  ISO  1=1.19 
150  01(11=02(1) 

200  CONTINUE 
RETURN 
END 
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CH 
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CH 
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CH 

CH 


PROGRAM  A AR( INPUT • OUTPUT  t PUNCH, TAPE5= 1 NPUT  » T APe6=OUTPUT . 
lTAPe7=PUNCH  »TAPei, TAPe2.TAPe3> 


THIS  PROGRAM  DECODES  METEOROLOGICAL  DATA  OBTAINED  «ITH  THE 

USA  AUTOMATIC  MET  SYSTEM  AT  THE  MHITE  SAND  MISSLE  RANGE  DURING  OCT 

/NOV  1974,  (REVISED  FOR  UNIVAC.JUNE  1975) 


20  MAY  75 


PROGRAM  AAR  READS  VARIOUS  INPUT  DATA  DESCRIBED  AS  COMMENTS  IN 
IN  PROGRAM  AND  COMPUTES  A FINAL  FORECAST  CURVE  ( HE IGHT , VI RTUAL 
TEMP  IN  OEG.C.  .U  AND  V COMP.)  FOR  A SPECIFIED  FIRING  TIME  AND 
FIRING  LOCATION 


DOUBLE  PRECISION  XX , Y V • XL, YL, RAO 
INTEGER  ABCOEFI 14) , GH I JKL ( 1 0 ) • MNOPQR 1 1 0 ) 

INTEGER  H1(7.7).T1(7. 7),Ul(7,7).Vt(7.7) 

DIMENSION  HH(49.4,2)  .TV(  49.4.2  ).U(  49.4*2  ) .VIAg*  4*2)  (XMILSK  21  ), 
1HHFNAL(4.2) .TVFNAL (4. 2) .UFINAL(4.2) .VFINALI 4. 2 ) * XNOTSl ( 2 I ) 
DIMENSION  JSAV (5)*TlM(S)*COISl2I) .XNOUC  2 1 ) * XNOV ( 21 1 • XNOTCf 2 1 ) . 
lXNOT(21 ).XNOH(21)*XNaTS2( 2 1 I * XMIL S2( 2 1 1 • 0 IF SEC( 21  I 
01 MENSION  UF(21 ) .VF (21 ) .TF(2t ). HF(21 )*PF(21 ) .PFOI S ) . X( 7 ) * V( 7) 
DIMENSION  UR( 12. 12) • VR( 12. 1 2 ) .HR ( 1 2 * 1 2 ) • TRI 12*12) 

DIMENSION  IUM( 147 ) . XS 4 25) . YS ( 25 ) . ES( 25 ) . 1 MO ( 12).PR( 12. 12) 
DIMENSION  HHTFGIS) .STV(S) . HHTF 4 5 ) .T VTF ( 5 ) .UTF ( 5 ) . VTF( 5 ) 

COMMON  /FMT/  XLAB 
EQUIVALENCE! I 7. NX ) . ( J7.NV) 


DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 


STV/1 1 1 .0, 1457.0.3012. 0.5574.0/ 

ID IF, I GMT/20, -420/ 

L1*L4,L12/1*4. 12/ 

I NO* YM AX* YM IN/21 *300.0. 1000,0/ 
JSAV/10.1.9.10.10/ 

TIM  /I  700.* 1830. .2000. .2130. ,2300./ 

PFO/870. 0*850.0*  700.0*500.0*0.0/ 

I MD/0, 31. 59* 90*1 20* 151. 181 •212.243*273*304*334/ 


FORMAT! IX *1 9H  NO  DATA  FOR  TF(REs*F7.1) 

FORMAT! IHl) 

F0RMAT(8F10.2) 

FORMAT! 1 0X.3F10.2) 

FORMAT  413A6.A2) 

FORMAT! /65H  LATITUDE  LONGITUDE 

1 */,2F15.2) 

FORMAT!/! 5X*9HNaRTHING=*D20. 8*1  an  EAST  1 NG= • D20 .8  I 

FORMAT! lX.I 2.F7.2*F10.2.F5.2*I5*12X*F8.2*F10.2,F5.2*(5) 
FORMAT !5X *215* 2F1 0.2 *2 110) 


PRINT  3 
PRINT  1100 
PRINT  3 


THE  NEXT  DATA  SET  IS  INPUT  READ  FROM  DATA  CAROS  (VARIABLES  ARE 
DESCRIBED  IN  COMMENT  CAROS  THAT  FOLLOW 
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RFAO(5.10)  ABCOEF 
«>R1TE(6.10)  ABCOEF 

REA0(5«30 I IQ*  XF 1 RE . XUL • XO . NX • YF IRE « YUL  t Y D . NY 
«RI  TEI  6*3  0)  IQ.XFIRE.X(JL.XO.NX*YFIRE.YUL*YO.NY 
PRINT  1002 
REAOIS.IO)  GHIJKL 
MRITE(6*10>  GHIJKL 

REA0(5*42 )NT*NL*TFCST*TFIRE* lOATE* JSTAT 
WRI TE (6*42INT*NL*TFCST,TFIRE*IOATE. JSTAT 
( PRINT  1002 

; PRINT  1002 

; WRI TE(6« lOOS) 

} PRINT  1002 

REAOIS.IO)  MNOPQR 
WRI TEI 6* 10)  MNOPQR 

RE AO  S.{XS(J).YS(J)*ESIJ)*J=1.JSTAT) 

PRINT  S.IXSI J)«YSCJ). ESI J).J=1* JSTAT) 

r 

f 

j CH  A8COEF=HEAOER  LABEL  CARO  FOR  FOLLOWING  DATA 

! CH  IQ=INDICATOR  SHOWING  HOW  FINAL  OATA  IS  TO  BE  PRINTEO 

' CH  1 PRINTS  HARO  COPY  OF  CURVES  ANO  OATA  FOR  CURVES 

CH  -1  PRINTS  CURVES  ONLY 

CH  XFIRE=LONGITUOE  OF  CENTER  OF  FIRING  PATH  (OEGREES  ANO  TENTHS) 

CH  XUL=LONGI TUOE  OF  UPPER  LEFT  HANO  CORNER  OF  GRID 

CH  X0=0ISTANCE  (OEGREES  ANO  TENTHS)  BETWEEN  GRIO  POINTS  IN  X (LONG.) 
CH  NX=NUM8ER  OF  MINI  GRIO  ARRAY  POINTS  IN  X OIRECTION  (LONGITUOE) 

CH  YF IHE=LAT ITUOE  OF  CENTER  OF  FIRING  PATH  (OEGREES  ANO  TENTHS) 

CH  YUL=LATITUOE  OF  UPPER  LEFT  HANO  CORNER  OF  GRIO 

CH  YO=OISTANCE  IN  OEG.  ANO  TENTHS  BETWEEN  GRIO  POINTS  IN  Y OIR.  (LAT) 

CH  NY=NUMBER  OF  MINI  GRIO  ARRAY  POINTS  IN  Y OIRECTION  (LATITUDE) 

CH  GHl JKL=HEAOER  LABEL  CARO  FOR  FOLLOWING  OATA 

CH  NT=NUMBER  OF  I HR  FCST  ARRAY  PERIODS  INCLUDING  TIME  OF  FCST 
CH  NL=NUMBER  OF  LEVELS 

CH  TFCST=GMT  TIME  OF  FIRST  ARRAY  OF  UPDATED  GWC  FORECAST 

CH  TFIRE=GMT  TIME  OF  FIRING  (HRS  AND  MIN) 

CH  IDATE=DATE  OF  GWC  FORECAST 

CH  JSTATsNUMBER  OF  OBSERVATION  STATIONS 

CH  XS( J)-LONGITUOE  OF  REPORTING  STATION  IN  OEG  AND  HNOTHS 

CH  YS( J )=LAT ITUOE  OF  REPORTING  STATION  IN  OEG  ANO  HNOTHS 
CH  ESI J)=ELEVATION  OF  REPORTING  STATION  IN  METERS 

IY=INT(  IOATE40.0001  ) 

IM=INT( IOATE40.01 )-IY*100 
ID=IDATE- IY4 10000-IM4 I 00 
IH=INT(TFIRE) 

IMN=( TFIRE-IH) *60 

IMOATEs  IMN'feO*!  IH4-24*(  104  I MO  ( I M ) 4 365*  I Y-366  ) ) 

REX=0.0 

nExi=o.o 

XEND=TFCST4NT-I 

CALL  TIME  (TFIRE.TFCST.TS.TL.ITS.ITL.XEND.REXl) 

IF  (REX1.E0.59)  GO  TO  600 

CALL  CENTER(XFIRE*YFIRE*XUL*YUL*XO*VD«NX*NY*PN*PE«PS«PW. 

1 IPN.IPE. IPS.IPW.REX) 
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r ^ 

[ IF  (REX.EQ.71)  GO  TO  600 

[ ITSS=ITS-I 

r 

I 

CH  THE  FCST  DATE  (JOATE)  ANO  TIME  IJTIME)  OF  G«C  DATA  IS  NOW  READ  IN 

REAOIll  JOATE. JTIME 

! XTIME=FLOATl JTIME)*O.OI 

[ IFIJOATE.NE. lOATEIPRINT  1001 

I IFIXTIME.NE.TFCST ) PRINT  1001 

I CHECK  IFIJOATE.NE.IOATEIGO  TO  600 

CHECK  IFIXTIME.NE.TFCST)  GO  TO  600 

IF  (TFIRE.GT.XENO.OR.TFIRE.LT.TFCST)  PRINT2.TFIRE  , 

CHECK  IF  ITFIRE.GT.XEND.OR.TFIRE.LT.TFCST)  GO  TO  600 

CH  THE  GWC  FCST  DATA  IS  SCANNED  AS  I DIM! ) TO  DETERMINE  THE  FCST  DATA 

CH  ON  EACH  SIDE  OF  THE  T I ME  OF  FIRING 

f 

PRINT  42.ITSS 

IF!  ITSS.EQ.O)  GO  TO  202 

DO  200  JT=1.ITSS  i 

DO  200  JL=L1*L4 
200  READ!  1 HUM 

CH  IN  THE  211.210  DO  COOPS.  THE  GWC  FCST  DATA  FOR  LEVELS  LI  TO  L4  IS 

CH  NOW  READ  FOR  THE  FCST  TIME  ON  EACH  SIDE  OF  THE  TIME  OF  FIRING 

202  DO  211  JT=1  .2 

DO  210  JL=L1.L4 

REAO(l)  ( lUllJ. 1 I .VI ( J.l) .HICJ. 1). I3l. 17>.Jsl. J7) 

DO  205  1=1.17 
00  205  J=1.J7 
1H1=1NT(H1IJ.I )«0.0001) 
r 1 fj. 1 )=1ABS( HI fj. I )-lM 1*1 0 000) 

HlfJ.l)sIHl 
205  CONTINUE 

DO  209  JX=1.I7 

DO  208  JV=l.J7 

JXV=JX4-(  JV-1  )*NX 

HH( JXV. JL. JT)=H1 (JX. JY)*0. 1 

TV(JXY.JL.JT)=Ti( JX.JY)*0. I 

U( JXY.JL.JT)=U1 I JX. JY)«0.1 

208  V( JXY. JL. JT)=V 1 I JX. JY) *0. 1 

209  CONTINUE 

210  CONTINUE 

211  CONTINUE 
JX=1PW 
JV=IPS 

JXY=  JX*-(  JY- 1 ) ♦NX 
JWS=JXY 

J X=  I PE  ■ ^ 

JXY  = JX  + I JY- 1 HNX 

JES=JXV 

JV=IPN 

JXY=JX»I JY-l )*NX 

JEN=JXY 

JX=IPW 


JXV= JX* ( JV- I ) *NX 
JWN= JXY 


CH  IN  THe  350  00  LOOPS.  HEIGHTS.  TENPERATUHES.  U AND  V COMPONENTS  FOR 

CM  THE  VARIOUS  LEVELS  IJL=LI.L4)  ARE  COMPUTED  FOR  THE  LOCATION 

CM  OF  XFIRE  AND  YFIRE 

DO  250  JT=I.2 
DO  250  JL=L1.L4 

CALL  INTERPIPM.PE.HHl  JVS. JL. JT)  .HH I JES . JL  . J T ) .XFIRE. 

IHMXFS) 

CALL  I NTERPIPW .PE.HHI JVN.JL.JTI.HHIJEN.JL.JT). XFIRE. 

1 HHXFN) 

CALL  I NTERPIPS.PN.HHXFS.HHXFN.YFIRE.HHFIN) 

HHFNAH  JL  . JT)  = HMFIN 

CALL  I NTERPIPM .PE.TVIJWS. JL. JT  ). T V( JES. JL • JT ) .XFIRE. 

ITVXFS) 

CALL  1 NTERPf PM. PE. TV  I JVN.JL .JT ).TV( JEN . JL • J T ) . XF IRE. 

1 TVXFNI 

CALL  INTERPIPS.PN.TVXFS.TVXFN.YFIRE.TVFIN) 

T VF  NAL I JL  . JT  ) =T  VF  I N 

CALL  INTERP(PW.PE.U( JWS.JL. JTI.UIJES. JL. JT) .XFIRE.UXFS) 

CALL  I NTERPf PW.PE.UI JWN. JL . JT  ).U( J EN. JL . JT ) . XF I RE . UXFN ) 

CALL  I NTE RPIPS.PN.UXFS. UXFN. YF IRE .UFIN) 

UF  INAL!  JL.  JT  »=UFIN 

CALL  I NTERPIPM.PE.V ( JW  S. JL. JT). VIJES.JL . JT ) .XFIRE. VXFS) 

CALL  I NTERPIPW.PE. V( JWN. JL. JT) . VI  JEN. JL. JT ) .XFIRE.VXFN) 

CALL  INTERPIPS.PN.VXFS.VXFN. YF IRE. VFIN) 

VF INAL! JL. JT)=VFIN 
250  CONTINUE 

CH  IN  THE  260  DO  LOOP.  HEIGHTS.  TEMPERATURES.  U AND  V COMPONENTS  FOR 

CH  THE  VARIOUS  LEVELS  (JL=LI.L4)  AT  LOCATION  XFIRE. YFIRE  ARE  COMPUTED 
CH  FOR  THE  TIME  OF  FIRING  (TFIRE) 


260 

435 


DO  260  JL=LI.L4 

CALL  INTERP  ITS .TL . HHFNAL I JL . 
CALL  INTERP  I T S. TL . T VF N AL I JL . 
CALL  INTERP  I T S .TL .UF I NAL I JL . 
CALL  INTERP  I TS . TL . VF I NALI JL . 
CONTINUE 

DO  435  LHT=L1.L4 

HHTFGILHT )=HHTF( LHT )FSTVILHT ) 
PRINT  3 

WRITE(6. 1045)  JDATE.JTIME 

PRINT  1110 

PRINT  1112. UTF 

PRINT  1112. VTF 

PRINT  1112. TVTF 

PRINT  1112.HHTFG 


1 ).  HHFNAL  I JL  .2) 
1 ).TVFNAL( JL.2) 
1 ).UFINAL ( JL.2) 
1 ) .VF INAL I JL.2 ) 


-ESI  I A) 


. TFIRE.HHTF  I JL)  ) 
.TFIRE.TVTFI JL) ) 
.TF  IRE.UTFI JL ) ) 
.TF IRE.VTFI JL) ) 


CH  SOUNDINGS  OF  THE  AVAILABLE  U.S.  ARMY  RAWINSONOE  STATIONS  ION  TAPE) 

CH  ARE  READ  IN  IN  SUBROUTINE  RAWIN 


1 = 1 

JT=JSTAT 
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CALL  RAWI N( 1 . JT . 01 FS EC • ES . UH« VR t HR.PR. TR* IM  C. IMOA TE t 10  I F , IGM T ) 

CH  IN  DO  LOOP  311*  THE  UEST  AVAILABLE  STATION  IS  LOCATED  BASED  ON  THE 

CH  TIME  AND  LOCATION  OF  THE  SOUNDINGS  AVAILABLE 

DO  311  J=1.JT 

XL00=(XS<  J»-XFIRE)*1I 1«I37*COS(VFIRE*0.01745) 

YLAD=<  VS( J)-VF IRE) «1 1 1 • 137 

CO  I S(  J ) = AaS(OIFSEC(  J)  M.5*SQRT(  XLaO«XLaO'»’YLAO*YLAO) 

311  CONTINUE 
IA«1 

CMIN^t .E^OG 
DO  320  J=1«JT 

IFICOISI J ) .GT .CMlN.OH.UR(9. J ).LE.0.0 ) GO  TO  320 
CMIN=C0IS(J) 

I A=J 

320  CONTINUE 

CH  AT  THIS  POINT.  THE  CURVE  OF  THE  BEST  AVAILABLE  SOUNDING  IS 

CH  PRINTED  AND  THE  CURVES  ARE  MAPPED  FOR  BOTH  THE  BEST  AVAILABLE  AND 

CH  THE  UPDATED  SOUNDING 

DO  400  L=1.L12 
UF<L)=UR(L. lA) 

VF( L)=VR(L. lA) 

TFID^TRIL.  lA) 

HF(L)=HR(L. lA ) 

PFIL)=PR(L. I A) 

400  CONTINUE 
X( 1 1=XUL 
Yl 1 )=YUL 
DO  410  1=2.7 
X(  I J = X( I-l )4XD 
YI  I ) = Y*  l-I  )-YD 
410  CONTINUE 
XMAX=20. 

XMlN=-20. 

CALL  MOVCURCUF.PF.UTF.PFO.XNOU) 

IF! IQ.LT. 1)  GO  TO  420 

XLAB=6HU-C0MP 

PRINT  1010 

CALL  PNTDATIUF.PF.XNOU.PF. IND. INO. XMAX . XM IN . YMA X . YM IN  . SOB  I 
420  CALL  MOVCURIVF.PF.VTF.PFD.XNOV) 

IF( IQ.LT. 1 ) GO  TO  430 

XLAB=6HV-C0MP 

PRINT  I 01  I 

CALL  PNTOATI VF.PF.XNOV.PF. IND. fNO. XM A X. XM IN . YMA X . YM I N . BOB ) 

430  CALL  MOVCURITF .PF.TVTF.PFO. XNOT ) 

IF( IQ.LT. 1 ) GO  TO  440 
XMAX=31 0. 

XM1N=21  0. 

XLAB=6H  TEMP 
PRINT  1012 

CALL  PNTOATI TF.PF.XNOT .PF . INO • I NO. XMAX . XM IN . YMAX . YM IN. BOB ) 

440  CALL  MOVCURIHF.PF.HHTFC.PFO.XNOH) 

IFIIQ.LT.I)  GO  TO  450 
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XMAX=I6000. 

XMI N=0. 0 
XUAa=6HHe  IGHT 
PRINT  1013 

CALL  PNTDAT (HF .PF.XNOH«PF.lNO.IND«XMAX«XMIN.YMAX.yMIN«BaB) 

450  CF(UF( t ) .FO.O.U .ANO.VFI I ) .EO.O.O ) GO  TO  460 
GO  TO  500 

460  IFITFI  1 ) .EQ.O.O.AND .HF  ( 1)  . EQ.O. 0)  PRINT  2.TFIRE 
500  CONTINUE 

C *STM  CONVERSION 

RAD=0.0 1 74532925 
YL=rF  IHE4RAO 
XL=XF IBE4RAD 
CALL  WSTMIYL.XL.XX.YY) 

PR  I NT  3 
PRINT  1046 
PRINT  1002 
PRINT  1040 
PRINT  1002 
PRINT  1020 
PRINT  1002 
PRINT  1025 
PRINT  1002 

CH  IN  00  LOOP  470.  DATA  FOR  BOTH  THE  BEST  AVAILABLE  SOUNDING  AND 

CH  UPDATED  SOUNDING  ARE  PRINTED  (F.l  FORMAT)  FOR  ZONES.  PRESSURE, 

CH  temperature,  and  U and  V COMPONENTS 

DO  470  NO=l . 12 

IF  (PF(NO).GT.O.OI  GO  TO  523 

TF(ND)=-999.9 

UF ( ND)=-999.9 

VF( ND) =-999.9 

XN0T(ND)=-999.9 

XNOU(NO)  =-999.9 

XNOV(NO)  =-999.9 

523  PRINT  1030.  ND ,PF ( NO ) . TF { NO ) . UF ( NO ) . VF ( NO) . ND . PF{ ND ) , 

1 XNDT ( NO) . XNOUI NO) .XNOVI ND ) 

470  CONTINUE 

PRINT  1002 

PRINT  12.  YFIRE.XFIRE 

PRINT  3 

PRINT  1050 

PRINT  1002 

PRINT  1060 

PRINT  1002 

PRINT  1020 

PRINT  1002 

PRINT  1070 

PRINT  1002 

PRINT  1002 

CH  IN  DO  LOOP  460.  DATA  FOR  BOTH  THE  BEST  AVAILABLE  SOUNDING  AND 

CH  UPDATED  SOUNDING  ARE  PRINTED  IN  COMPUTER  MET  MESSAGE  FORMAT 
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475 


4 76 
4 78 


533 

480 


tool 
1002 
1005 
1010 
1011 
1 01  2 
t 01  1 
I an 
I 025 

I 1 10 

• 1 «l 
' *«  I 

♦ •4 

»•  » 


DO  480  N0=1 ,12 

IFIUFINDI  .EQ.O.O.ANO. VF(NO) .EQ.O.OI  GO  TO  4 74 
DIRl  =ATAN2(-*JF1N0  ) .-VF  ( NO}  ) 4 180/3.  14  159 
IF  lOIR 1 .LT .0.0 ) OIR13DIR14360 
XMILSl ( NO }s01R 14(6400. 0/360  .0  >/10.0 
XMPSl«SaRT(UF(N0)4UF( N0)4VF( NO ) 4VF ( NO > ) 

XNOTSl (N0I«1.942S44XMPS1 
GO  TO  475 
XMILSKNO)  a-999 . 9 
XNOTSl (N0)«-999.9 

IF( XNOU(NO) .EQ.O. 0. ANO.XNOVf NO) .EQ.0.0)  GO  TO  476 

0IR2:::AT  AN2(-XN0U(N0  ).-XNOV(ND)  14180/3.14159 

IF  (OIR2.ET.O.O)  OIR220IR2'<-360 

XM ILS2( NO )301R24( 6400. 0/360.0)/ 10. 0 

XMPS2  = SQRT(XNOU(NO)4XNOU(NO)4-XNOV(NO)4XNOV(NO)  ) 

XNOTS2(  NO  )-=  1. 942544  XMPS2 

GO  TO  478 

CONTINUE 


N=NO-l 


IXNOT=IFIX( 10. 04X NOT (NO >4. 5) 

IPF=lFIX(PF(NO)4.S> 

I XMIL1=IF(X (XMILSl I NO >4.5 ) 

I XNOTl  = IFI X( XNOTSl I NO  >4.5) 

IXMIL2=1F1X(XM1LS2(N0)4.5) 

I XNOT2=IF IX( XN0TS2( N0)4.5 ) 

ITF^IFIXI 10.04TF( NO)40.5) 

IF  (PF(NO).GT.O.O)  GO  TO  533 

ITF  =-999 

IXN0T=-999 

I XN0Tl=-999 

IXNOT2=-999 

IXMILl=-999 

I XMIL2-=-999 

PRINT  lOaO.N. IPF. ITF.  IXMILl . IXNOTl.N.  IPF.  IXNOT. IXM1L2. I XN0T2 

WRITE  (3)  N.  IPF  . I TF  ,1  XMILI  • IXNOT  1 *N.  IPF.  IXNOT.  IXM1L2.  1XN0T2 

CONTINUE 

PRINT  13.YY.XX 

WRITE  (3)  XX. YY 


OATA) 

CF  GRAPHICAL 
OF  GRAPHICAL 
CF  GRAPHICAL 
OF  GRAPHICAL 


OATA 

OATA 

OATA 

OATA 


FORMAT! 1X.46H0ATE  TIME  OF  GWC  FCST  OlFFERENT  THAN  OATA  CARO) 
FORMAT! IX. 21 15) 

FORMAT ( 1 7X. 17HARMY  STATION 
FURMAT!  IHl  . 15X  .46HPRI NT CUT 
FORMAT! IHl . 15X.46HPR(NT0UT 
FORMAT! 1 HI . 15X .44HPR1 NT CUT 
FORMAT!  IHl . 15X .46HPRI NTOUT 
FORMAT (4X.23HeEST  AVAILABLE  SOUNO I NG. 20X. 16HUP0 ATEO  SOUNO(NG) 
FORMAT!! X.70HZ0NE  PRESS  TEMP  UCOMP  VCOMP  ZONE  PRESS 

IMP  UCOMP  VCOMP) 

F0HM4T( 2X.1 3.F7.1 .FG. 1 .2F7. 1 .1 1 1.F7. 1 .F6. 1.2F7. 1 ) 

FUMM4TI // 10X.5HCM1N=.F 10.2. lOX. 3Hi A=. IS) 

FOMMwrt 7K.49HPRESS  IN  MBS.  TEMP  IN  OEGlK).  U ANO  V COMP  IN  MPS > 
F«>mm4TI  lOX.  lOHGwC  OATE  =.  I 1 0 . 1 OX.  1 OHGWC  TIME  =.110//) 

F MMaAT I I6X. ISMSTANOARO  FORMAT) 

• •••r 12 1 I.ZZMCOMPUTER  MET  MESSAGE  FORMAT) 


(X=U  COMP.  Y=PRESS)/) 
!X=V  COMP.  Y=ORESS)/) 
(X*TEMP.  Y=PRESS)/) 
IX=HEIGHT.  Y=PRESS)/) 


TE 
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1060  FORMAT ( IX .8  IMPRESS  IN  MHS t TEMP  IN  TENTHS  OF  OEGlKl*  DIRECTION  IN 
ITENS  OF  MILS.  SPEED  IN  KNOTS) 

1070  FORMAT! IX. 70HZ0NE  PRESS  TEMP  DIR  SPEED  ZONE  PRESS  TE 

IMP  OIR  SPEED  ) 

1080  FORM  AT  I 2X. 13. 17. I6.217.I11. 17. 16.217) 

1100  FORMAT! ///ATX, 30MARTILLERY  APPLICATIONS  ROUTINE.  ////// 

2///.  57X.5HUNITS.///.53X. I IHSPEED  - MPS. // . 53 X. 1 5HO I REC T I ON  - OE 

3G.//.  53X. 15HHE1GHT  - METERS. //. 5 3X . 1 OHTEMPERATURE  - DEG  K) 

1110  FORMAT (bX.SHGMC  DATA//) 

1112  FORMAT (4X. lOEl 2 .3/ ) 

600  STOP 
END 


84 


SUBROUTINE  TIME  I T F 1 HE  * TFC ST t T 5 . TL . I TS . 1 TL . TFENO . RE X I ) 

CH  SUeWOUTINECTIME)  SELECTS  THE  MINI  GRID  FORECAST  ON  EACH  SIDE  OF 

CH  THE  TIME  OF  FIRING  AND  INTERPOLATES  TO  TIME  OF  FIRING 


CH 

CH 

CH 

CH 

CH 

CH 


TFIRE=TIME  OF  F 
TFCST=  TIME 

TS=TIME  OF  FCST 
TL=TIME  OF  FCST 
ITS=ARHAY  POSIT 
ITL=ARRAY  POSIT 


IRING 


I HRS  AND  DECIMALS) 


OF  FIRST  ARRAY  OF  UPDATED  GI«C  FORCAST 

ARRAY  THAT  IS  JUST  SMALLER  THAN  TFIRf 

ARRAY  THAT  IS  JUST  LARGER  THAN  TFIHE 

ION  INTEGER  OF  TS 

ION  INTEGER  OF  TL 


20  F0HMAT(5X.37HTFIRE  IS  GREATER  THAN  FORECAST  EXTENT) 
IF  ( TF  IRE.LT.TFCST)  TF  IRE=  TF  IRE4-2 A . 0 
IFITFIRE.GT.TFENO)  REX1=S9 
IF( TF IRE, GT. TFENO)  GO  TO  40 
ITS  = 1 
TS=TFCST 
DO  80  ITL=2*19 
TL=TFCSTaITL-1 

IF(TFIRE.LT.TL)  GO  TO  100 


I TS=ITL 

80 

TS=TL 

REXl=59 

40 

WRI TEIb 

100 

RETURN 

END 
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SUBROUT  I NR  CENTER  ( XF  I R E t VF  IR E , XUL  t YUL  . X D » VO . NX » NY , PN . PE  » PS,  PVK, 
IIPN, IPE, IPS.IPK.REX) 

CH  SUBROUTINE  CeNTEH=SCANS  JX  AND  JV  ARRAYS  AND  INTERPOLATES  VALUES 

CH  FOR  THE  POSITION  XFIRF.  AND  YFIRE 


CH  XF IHE=LONGITUOE  OF  CENTER  OF  FIRING  PATH  (DEGREES  AND  TENTHS J 

CH  YF IRE=LAT ITUOE  OF  CENTER  OF  FIRING  PATH  (DEGREES  AND  TENTHS) 

CH  XUL=LONGI TUDE  OF  UPPER  LEFT  HAND  CORNER  OF  MINI  (WIO  ARRAY 

CH  YUL=LATITUDE  OF  UPPER  LEFT  HAND  CORNER  OF  MINI  GRID  ARRAY 

CH  X0=DISTANCE  (DEGREES  AND  TENTHS)  BETWEEN  GRID  POINTS  IN  X OIHECTION 

CH  YO=DISTANCE  (DEGREES  AND  TENTHS)  BETWEEN  GRID  POINTS  IN  V DIRECTION 

CH  NX=NUMBER  OF  MINI  GRID  POINTS  IN  X OIHECTION  (LONGITUDE) 

CH  NY=NUMBER  OF  MINI  GRID  POINTS  IN  V DIRECTION  (LATITUDE) 

CH  PN=POSITION  (LATITUDE)  OF  THE  ARRAY  GRID  JUST  NORTH  OF  YFIRE 

CH  PE=POSITION  (LONGITUDE)  OF  THE  ARRAY  GRID  JUST  EAST  OF  XFIRE 

CH  PS=POSITION  (LATITUDE)  OF  THE  ARRAY  GRID  JUST  SOUTH  OF  YFIRE 

CH  PW=POSIT10N  (LONGITUDE)  OF  THE  ARRAY  GRID  JUST  WEST  OF  XFIRE 

CH  IPN- ARRAY  GRID  NUMBER  JUST  NORTH  OF  YFIRE 

CH  IPE=ARHAV  GRID  NUMBER  JUST  EAST  OF  XFIRE 

CH  IPS- ARRAY  GRID  NUMBER  JUST  SOUTH  OF  YFIRE 

CH  IPW=ARRAY  GRID  NUMBER  JUST  WEST  OF  XFIRE 


100  FORMAT  (1X,29H 

110  FORMAT  (1X,29H 

120  FORMAT! IX, 2BH 
130  FORMAT! 1 X ,2aH 


YFIRE  NORTH  OF  MINI  GRID) 
YFIRE  SOUTH  OF  MINI  GRID) 
XFIRE  WEST  OF  MINI  GRID) 
XFIRE  east  of  mini  GRID) 


PE=0.0 
PW=0.  0 
PS=0.0 
PN=0,0 

IFIYUL.LT.VFIRE  ) PRINT  ICO 
YLL=YUL-YO*NY 

IF( YLL.GT. YFIRE ) PRINT  110 
IFIXUL.GT. XFIRE)  PRINT  120 
XUR=XUL+XD*NX 

IF ( XUR.LT ,XF  IRE  ) PRINT  130 
DO  200  I X=1  ,7 
PE=XUL*XDF ( I X-1 ) 

IF  (PE.LT.XF IRE)  GO  TO  200 

1PE=  IX 

pw=PE-xn 

IPW=IPE-1 

GO  TO  201 

200  CONTINUE 

201  DO  220  IY=1 ,7 
PS=YUL-VO*( lY-l ) 

IF(PS-GT.VF IRE ) GO  TO  220 
IPS=I Y 

PN=PSFYD 
I PN=I PS-I 
GO  TO  222 
220  CONTINUE 
222  RETURN 
END 
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SUBROUTINE  INTERPI X 1 • X2* Y 1 • V2 . XO . YO ) 


4 


CM 

CM 


SUBROUTINE  INTERP  INTERPOLATES  INTERMEDIATE  VALUE  OF  YO  FROM  INITIAL 

XI.  VI .1 NTERMEOI ATE  XO  AND  FINAL  X2 » Y2  VALUES 


I F(  Y2.LE.VI  »Y0=<  I X2-X0  )*I  YI-V2  »/<  X2-X  I ) )4-Y2 
IF  I VI  .LT.Y2>Y0=( IXO-Xl  ) ♦< Y2-Y I ) / ( X2-X 1 » )+Yl 
RETURN 
END 


nnnnnr*  nr> 


SUOROUTINE  RAWiNd  • JT. 7S«ES« UR • VR.HR.PR.tr*! MO t I OATE.IDIF,  IGMT ) 
THIS  SUbKOUriNE  READS  IN  THE  U.S.ARMY  RAWINSONOE  DATA. 

DIMENSION  UR(  1 ) . VR ( 1 ) .HR ( 1 ) .TR( t ) .PR( 1 ) .ES( 1 ) .TS( 1 ) . IMO(  I ) 
DIMENSION  O AT  I 512 ) • IS (25 ). HOZ( |7  > 

DATA  ACR.CKM.XNlL/0.0 174533. 0.5 I 4n.-999.9/ 

DATA  J3.J10.J12.J17.J23.J28.J5 1/3 .10. 12. 144. 23. 28. 51/ 

DATA  HOZ/0.0. 1 00* 0*350. 0.750.0 . 1 250- 0 . 1 750 • 0. 2250. 0. 2 750 .0 . 3250 . 0 
2 . 3750 .0. 4250.0.4750.0.5500. 0.6500.0. 7500. 0.8500.0.9500.0/ 

I MI N= I 44  0 
DO  18  J=1.JT 

18  1S(J)=0 
DO  19  J=1 . J1 7 
HR<  J)=XNIU 
TR( J)=XNIL 
PHI J)=XNIL 
UR( J)=XNIL 

19  VR(J)=XNIL 
IF  ( I .GT  .1  » GO  TO  22 

HEADS  IN  DATA  SET 

20  READ  2. DAT 
2 FORMAT ( 1 OFS. I ) 

LOOP  THROUGH  DATA  SET  IN  RECORD 

21  DU  30  J=1.J10 

22  JC=(  J-l)*J5l-H 
C 

C CHECKS  DATA  FOR  DA TE . T I ME . AND  TYPE 
C 

IP1=INT(0AT( jC)*1.0E-4  ) 

DAT J=DAT( JC)-IP1*I .0E4 
IP2=INT<DATJ41 . OE-2 ) 

OAT J = DATJ- 1P24 1 .0E2 
IP3=INT(0ATJ) 

IP4=INT(0AT( JC4I)*0.01 > 

I P5=DAT<  JC4  1 )-IP44100 

JOATE=  IP546  0*(  IP4424*(  |P34-1MD(  I P2  ) 4365*  1 PI -366  ) )-IGMT 
IF  ( JDATE.LT.IDATE-IMIN)  GO  TO  32 
IF  ( JDATE.lt . IDATE-IDl F ) GO  TO  30 
IF  I JDATE.GT. I0ATE410IF)  GO  TO  32 
24  JS=INT(OAT( JC42)*0.1  ) 

DAT2=OAT( JC42)-JS*1 0 

IF  (0AT2.LT.1.0R.0AT2.GT.3  ) GO  TO  30 
ISC  JSI  = INT<0AT2) 

IF  ( (S( JS  I .EQ.  1 1 GO  TO  30- 
ESJS=eS<  JS) 

TS<  JS)=JDATE-IOATE 

C SELECTS  OUT  PRESSURE. TEMPERATURES  AND  HINDS  FOR  ARTILLERY  ZONES 
JC  J- JC  + J28-2 
JKT=1 

DO  28  JH=l . J23.2 


86 


27 


28 

30 

32 


JK=(  JS-I  ) *JI2«'JKT 

JJH= JC JM 

HR<  JK )sHOZ<  JKT  I 

IFC ABS(DAT( JJH) ).GT.360.0I  GO  TO  27 
ANG=ACR«OAT( JJH) 

UR( JK)=-DAT(JJH* 

VR<JK)=-OAT<JJH«- 

IFCOAT( J3*JH).LT.0.0)  GO  TO  28 
PRi JK)sDATC  J3>JH) 

TR  ( JK  )=OAT(  J3-( 

JKT=JKT*I 
CONTINUE 
GO  TO  20 
CONTINUE 
RETURN 
END 


1 

) 

*S1N 

( AN6> 

1 

) 

*COS 

( ANG) 

• 

0 

.0  ) 

GO 

TO 

) 

4* 

1 

>♦27 

3. 

16 

* 

SUBROUTINE  MOVCURIXO. YC*XFD« VFO.XNO ) 

DIMENSION  XOI  I)»YO(  1 ) .XFO ( I ) • YFO ( I I . X NO ( 11 

DIMENSION  XfS) (EIS) • IY(5) 

C OETERMINE  THE  NUMBER  OF  FQFCAST  DATA  POINTS 
NPTSF=0 

DC  120  1=1*5  ^ 

120  IF! XFO( I ) .GT.O. 0)  NPTSF=NP T SF* 1 ' 

C DETERMINE  THE  NUMBER  OF  BALLOON  DATA  POINTS 
NPTSO=0 
DO  130  1=1.21 

IF  I Y01 I ) .LE.O.O.OR.XDin .LT.-900«0}  GO  TO  132 
130  NPTSD=NPTSD*1 
132  NPTSO=NPTSD-l 

C CALCULATE  X ORDINATE*  ERROR*  AND  INDEX  OF  BALLOON  DATA  RELATIVE  TO 
C FORCAST  DATA 

DO  no  J=I*NPTSF 
DO  140  I=t.NPTSO 

IF{  YD(  I l.GE.YFDI  J).AND-YD(I  + l).LT*YFO(jn  GO  TO  100 
140  CONTINUE 

IF  ( YFOIJ I .LT.YUI 1 ) ) GG  TO  142 
XIJ)=XD(NPTS0+1  ) 

IY( J)=NPTS0*1 
E( J)=XFO( J»-X< J) 

GO  TO  110 
142  X(J)  = XO(n 
lYC  J ) = 0 

E<  J)=XFD(  J)-X<  J) 

GO  TO  no 
100  CONTINUE 

X(  J ) = ( YD(  1 >-YFO<  J)  )*(  xon^l  )-XO<  I ) ) /(  YDl  I >- YD ( I ♦ 1 ) n-XD(  I ) 

E( J)=XFOI J )-X ( J ) 

IY( J)=I 
110  CONTINUE 
IS=IY( 1 ) 

IF  (IS.LT.l)  GO  TO  150 
DO  190  1=1* IS 
190  XNDI I ) = XO(I  )4E(  n 
150  NPTSF=NPTSF-I 

DO  200  J=1*NPTSF 

EF= (E( J4l »-E( J ) 1/1 YFDI J41 >- YFOIJ) ) 

IS=IY(jn-l 
DO  180  I=IS*NPTSD 

XND(1)=X0<I)*(YD(I )-YFO(J ) ) AEF+EI J ) 

IF! I .GE*I Y| J+1 ) )GO  TO  200 
180  CONTINUE 
200  CONTINUE 

IS= I Y( NPTSF+ 1 ) 

NPTS0=NPTSD4I 
DO  170  I=IS*NPTSO 
170  XNDI  I ) = XD(  I )4E(NPTSFF 1 ) 

RETURN 

END 

^Thls  subroutine  was  written  by  Len  Gaalorek  (SRI) 
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SUHftUUT  I NF  PNTDAT(xn,V0»Xi?tY2«N.M,XMAX,X^<IN.YMAX,Y'<II  N.TllH  \ 


■k 


c THl<;  PHOGRAV  PRCOUCES  NSFTS  NUMBER  OF  PR  I NT  F R-G  WA^HS  O-  X,Y 

c OAT  A RUUNneO  TO  THE  NF.AREST  OFL  X ANO  OELY. 


t I NPUT  . . . 

C XD  - X-OATA  FOR  FIR‘jT  SET. 

C YO  = Y-UATA  FOR  FIRST  SET. 

C XZ  = X-OATA  for  second  SET. 

c YZ  = y-oata  for  SEFUNO  set. 

c XMA  X , X M I N,  AX  . YM  I N = MAX  AND  MIN  OF  DATA. 

c N,M  = totals  of  o and  z type  data. 

C XLAlt.YLAO  = X AND  Y LABELS. 

C 

DIMENSION  XDl  l).YO(  I),X/(  l).YZ(  II 
COMMON  /FMT/  XLAH 

p format  < I H I 

4 format  ( I X . 3HX  =.I0EI2..T) 

5 FORMAT ( 1 X.3HY  =,10E12.3) 

<S  FORMAT  ( 4 X,  1 OEl  2 . 3 ) 

7 FORM  AT ( NX. I 9HBEST  AVAILABLE  DATA) 
a FORMaTINX, INHGWC  UPDATED  DATA! 


PRI  NT 

PR  I NT 

;> 

PRINT 

7 

PPI  NT 

c 

30 

PRI  NT 

4,  1 XD(L)  .L=l . 

10  I 

PR  I NT 

N,( XDILI ,L=1 1 

«N> 

40 

PR  I NT 

? 

PRINT 

5 . I YD( L ) .L= I . 

I 0 I 

PRINT 

6 . I YO(L) .L=l I 

. M 

PRI  NT 

P- 

PRI  NT 

n 

PRINT 

P 

bO 

PRINT 

4 . ( XZ (L  1 .L=l . 

10  I 

PR  I NT 

6 , (XZ(L)  ,L=1  I 

, Ml 

70 

PR  I NT 

2 

PRI  NT 

5,  (YZ(L),L=1  . 

10  I 

PRINT 

6. < YZ(L)  .L  = l I 

,M) 

IF  ( XMAX.NE .-0 .0  I GO  TO  1?0 
XMAX=-1 .0E300 
YMAX=- 1 . 0E300 
XMI N= 1 . 0F300 

YMi N= I . oeooo 

DO  100  I=1.N 

XMAX  = AMAXI  ( XMAX  «XU ( I ) I 

yMAX=AMAXl(YMAX.YD< I )) 

XMI N= AM  I Nl  ( XM  IN  ,X  D(  I ) I 
YMIN=AMI N1 ( YMIN.YDI  I ) ) 

100  CONTINUE 


BEST  AVAIliBLE  COPY 


* 

This  subroutine  was  written  by  Hlsao  Shigelshi  (SRI) 
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XMAX=AMAXl(XMAXtH/(l»)  , 

tYMAX  = A^'AXllY'«AX»Y'2(t)> 

XyiN=AMIM(XMIM.XZ(l»>  , 

YMIN=AMINl(YMIN.YZ( f)) 

I 1 n CONT  UJUf  ■ , 

1?0  CALL  PNT  M AP  (X  0 , YD  ,X  / , Y Z t Nt  M . XMAX  t XM  I Nt  AX  , YM  IN  ) 
j-  Pl'TUMN  I 

I f-ND 


! 

i 

t 

i 


n r 


1 

4 

*S 

ft 

c s 

1 0 

c 

c 

100 

c 


1 1 0 


I 20 

c 

1 004 


Sljn-?(IUTr  NF  PNTMAP(  AU.  VD.X/,  Y7  ,N,  M.XMAX  .XWI.N,  VMAX,  YM  IN  )* 

OIWtN'i  ION  XD(  1).Y0(  l).XZ(  l),Y7<  1) 

D I KIN  X ( 1 ) t Y(  ^ I ) , A(  ? 1 , Z 1 ) 

/F^IT/  XLAo 
DATA  YLA:?/ftH  PPESD/ 

FilPv^AT  (/IH  « 1 OX  , 1 3MX-  INCPt  Mt  NT  = . F 1 0 . T , X , 1 A H AND  Y- I NC -TF.  M F N T 
1 F 1 O . 3// 1 H . A6 ) 

FORMAT  ( IH  . fcPXAh) 

FncvATI/lX.Pft.l,  IX,  1H(  ,21  (2XA1  ) ) 

FORMAT! /IH  ,10X,1H*  ,?0(2X1H* )) 

FORMAT!  IH  ,F1  1 . 1 , AF  10.  1 ) 

FORMAT  I5X  , /‘SHI  N 1 HF  FOLLOWING  GRAPH  + = HtST  AVAILABLE:  SOUNDING  AN 
ID  X = GwC  UPDATtrO  SOUNDING) 

F CRM  AT  ! 1 H ) 


no  10  J = 1 , 2 1 
ARRAYS 

no  10  1=1,21 

A!  J,  I )=IH 

Obtain  xrangf,  yrange,  oelx  and  dfly. 

XR  ANGE=XMAX-XM  IN 
[^ELX  = XRANGE/  ! N-1  ) 

YRANr.L  = YMAX-YMI  N 
nELY=Y RANGE/ !N-1 ) 

DETERMINE  THE  X AND  Y VALUES  FOR  THE  AXES. 

X ! 1 ) = X M I N 

Y(  1 ) =Y  MI  N 

fVJ  100  I = 2,N 

X ! 1 )=X ! I -I ) FOELX 

Y! I )=Y! l-l)+OELY 

CONTI NUE 

DETERMINE  THE  IX,  lY  INDICES  FOR  ARRAY  A WHICH  HAS  THE  ITH  CO 
RELATED  TO  DELX  AND  JTH  ROW  RELATED  TO  UELY.  FOR  DATA 


SET  1 AND  FOR  SET  2,  INSERT 

AN 

X. 

DO  110  I =l ,N 

1 X= ! X0(  I 1 -XM IN  )/DELX  ♦ 

1 .5 

JX= ( YD ! I ) -YMI N I /DELY  ♦ 

1 .9 

IF!  1X.GT.21  .OR.IX.LT.  1) 

GO 

TO 

1 10 

IF!JX,GT.21.nR.JX,LT.l) 

GO 

TO 

1 10 

A!  I X, JX  )= IHF 

CONTI NUE 

DO  120  I = 1 , M 

IX= ( X2 ( I )-XM IN )/DELX  + 

1 ,S 

Ja=! YZ! 1 )-YMIN)/OELY  ♦ 

1 .5 

I F(  JX  .C.T  ,2  1 .OR  . JX  ,LT  . 1 1 

GO 

TO 

1 20 

IF!IX.GT.21.0H.IX.LT.l) 

GO 

TO 

120 

A ( I X, JX ) = IHX 
CO,\T  I NUE 

PRODUCE  THE  PRINTER  GRAPH  WITH  AXES  X AND  Y AND  ARRAY  A. 
PR  INT  1 OOA 
FORMAT ! I H 1 ) 

PRINT  S 


This  subroutine  was  written  by  Hlsao  Shigeishi 


(SRI) 
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PRINT  a 

PRINT  8 
PRINT  6 

PRINT  1 .OF-LX.OELY.YLAB 
I I =N^l 

DO  l-IO  I=I.N 
I 1-1  1-1 

PRINT  J.Y(II>t(A(J.IIJ.J 
CPNT I NUF 
PRINT  A 

PRINT  5t ( X( L ) .L=I . N .5 > 
PRINT  O.XLAB 
Rfc'  TURN 
END 


1 . N) 


BEST  AVAIUBIE  COPY 


SUHHaUTINK  *lSTM  {YLAT.XLON.  XX,  YV)* 

CUNve^^SlQN  FRLIGKAM  CQNVtHTS  Gt^PDHTIC  POGITICN  TO  WSTM 
YLAT  = NORTH  LAItTUOE  IN  RAO  I ANS 

XLON  = kkest  longitude  in  radians 

XX  = X (WSTM  EASTING)  - FT 
YY  = Y ( ASTV  NORTHING)  - FT 

|■)OUBLt  (JRF.CISION  XX  , YY  ,YLAT  , XLON,  NN,  AA,  OH,  N,  A,  H.  C,  D,  , V, 
I tSO.  LPGO,  KO,  SINEl,  DLAM,  II,  YrTERS.  P,  A6,  IV,  VV,  r S . 

? Ill,  KK.  CM,  tt,r)TANX 

data  AA,  HB  , SINEl,  MEIERS  / 6 3 7 tJ  i?0  6 . 400  , 6 3S6S  83  . SO  0 , 

I 0,48481  3ftHl  10O'i36D-b,  3 . .?H  0 S ^ 3 1 3 33  3 3 3 30  0 / 


CM  = 3e?eoc.oDO  * sinei 

DLAM  = CM  - XLON 
OLAM  = OLAM  / SINEl 
NK  = 3. ?820641 3800 

1 N = (AA  - BE)  / (AA  + PB) 

rsa  = (AA  **  ;?  - Bb  ♦*  ?)  / aa  ♦*  2 

2 FPSQ  = (AA  *•  2 - OB  *♦  2)  / BB  **  2 

A = AA  * (1,000  - N + 1,2500  ♦ (N  •*  2 - N *♦  3)  I 1,26562500  ♦ 

1 (N**4-N**5)) 

P = 1,500  • AA  * (N  - N **  2 + 0,87500  ♦ (N  *•  3 - N **  4) 

2 0,85R375D0  * N **  5) 

C = 0,P375D0  * AA  ♦ ( N ♦♦  2 ~ N ♦ A 3 ♦ 0,75D0  * ( N »*  4 - N * * 5 > > 

D = (35,000  / 48,000)  ♦ A A * ( N ♦ * 3 - N * * 4 3 0,687500  * N ♦*  5) 


C = 0,P375D0  * AA  ♦ ( N ♦♦  2 ~ N ♦ A .3  ♦ 0,75D0  * ( N »*  4 - N * * 5 > 1 
D = (35,000  / 48,000)  ♦ A A * ( N ♦ * 3 - N * * 4 3 0,687500  * N ♦*  S: 
E = 0,61523417500  * AA  * (N  **  4 - N ♦*  5) 

KO  = 0,06^600 

V = AA  / (1.000  - ESO  ♦ (OSIN(YLAT))  *♦  2)  0,5 

NN  = A • YLAT  - B * OSIN(2.0DO  • YLAT)  C * OSIN  (4,000  * 

I YLAT)  - 0 * OSIN  (6,000  ♦ YLAT)  EE*  OSIN  (8,000  * YLAT) 

NN  = NN  * KO 
62  P = OLAM  * 1 ,00-4 

II  = ((V  * DSIN(YLAT)  * OCOS  (YLAT)  ♦ SINEl  **  2)  / 2,000)  • KO 

1 *1,008 

DT ANX=nS I N( YLAT )/DCOS ( YLAT ) 

III  = ((V  * SINEl  **  4 *,OSIN  (YLAT)  * (OCOS  (YLAT))  **  3)  / 

1 24,000)  * (5,000  - (OTANX  ) •*  2 ♦ <J,000  * FPSO  ♦ (OCOS 

2 (YLAT))  **  2 ♦ 4,000  * EPSO  **  2 ♦ (OCHS  (YLAT))  **  4)  * 

3 KO  * 1,0016 

A6  = (P  **  6 * (V  * SINEl  ♦*  6 * OSIN  (YLAT)  * (OCOS  (YLAT))  ** 

1 5)  / 720.000)  * (61.000  - 58,000  * (OTANX  ) **  2 E ( 

2 OTANX  ) A*  4 E 270.000  * EPSQ  * I OCOS  (YLAT))  A*  2 - 

3 330.000  * EPSO  * (OSIN  (YLAT))  A*  2 ) * KO  A 1,0024 

85  IV  = V A OCOS  (YLAT)  A SINEl  * KO  * 1,034 

VV  - ((V  * SINEl  **  3 * (OCOS  (YLAT))  AA  3)  / 6.000)  * (1,000  - 

1 (OTANX  ) **  2 E EPSQ  * (OCOS  (YLAT))  **  2 ) * KO  A 

2 1,0012 

S5  = dabs  (P**5)  * ((SINE1**5  * V * (OCOS  (YLAT))  AA  5)  / 120.000 

1 A (6,000  - 18.000  * (OTANX  ) **  2 E (OTANX  ) AA  4 

2 E 14,000  A EPSQ  * (OCOS  (YLAT))  AA  2 - 58.000  * EPSO  * (OSIN 

3 (YLAT))  AA  ?))  * KO  * 1,0020 

ic 

This  subroutine  was  provided  by  ASL  of  the  U.S.  Army  Electronics  Command 


BBl  AVAILABLE  COPY 
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AO 

A0496d5 

B 

1 

1 

1 

H 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

END 

DATE 

EtLMCO 

3-78 

OOC 

Appendix  E 


PDRR  DATA  INPUT* 


This  Appendix  and  the  following  Appendices  (F,G,  and  H)  provide  an 
example  of  how  both  GWC  prognostic  information  and  Army  Meteorological 
observations  are  used  to  obtain  a Computer  Met  Message  for  an  artillery 
firing  for  a requested  time  and  location. 
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CAOOS  O 


CARD  E 


CAPOS  F 


BEST  AVAIM3LE  COPY 


INOUT  OATA  FOR  PORR  RniJTlNF 


7 - 1 O^. 70 


MONTH 
I 1 


STATION  INFORMATION  (Locations  of  Artillery  Met  Sections)’’ 


NUMBER 

LONGITUTE 

L AT  I TUTF 

FLFVAT ION 

1 

- 1 06*  32 

32*40 

1230.00 

(LC-36,  TSX) 

2 

- 1 06  • IS 

32 . 40 

1311  *00 

(Orogrande,  ORO) 

3 

- 1 06*00 

32*  27 

1 3S0. 00 

(Las  Cruces,  LSX) 

4 

- 1 06* 1 A 

32.27 

1 24Q.OO 

(MacGregor,  MCG) 

5 

-106*40 

32*27 

1 220 .00 

(War  Road,  WAR) 

6 

- 1 06*  48 

32*40 

1247.00 

(MTTR,  MTR) 

7 

- 1 06*  4? 

32*47 

1 220 .00 

(Small  Missile  Range,  SMR) 

9 

- 1 06*  IS 

32*S0 

1 2.30.00 

(Rampart,  RAM) 

0 

- I 06*  38 

32*62 

1210*00 

(Apache,  APA 

1 0 

-106*08 

32  *8S 

1250.00 

(Holloman,  HMS) 

♦Explanations  in  parenthesis  do  not  appear  on  actual  computer  output 
(See  Table  2 for  definitions  of  input  parameters) 
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- /^CPT 


Appendix  F 
AAR  DATA  INPUT 


I 

, 

[■ 

I. 


NT  NL  TFCST  IFIRE  lOATE  JSTAT 

19  4 0.00  12.30  741119  10 


ARMY  STATIUN  DATA 


LONG I TUDE 

LAT ITUOe 

EUEVAT I CN 

- 1 06. 32 

32.40 

1230.00 

- 106.  15 

32.40 

1311.00 

- 106. 90 

32.27 

1350.00 

- 106. 18 

32.27 

1249. 00 

-1 06.40 

32.27 

1220.00 

- 106. 48 

32.40 

1247.00 

-106.42 

32  .47 

1220. 00 

-106.15 

32.50 

1230.00 

- 106. 38 

32.  62 

1210.00 

-106. 08 

32.85 

1250.00 

BtSl'AVAIlABlE  COPY 


See  Table  3 for  explanation  of  Input  parameters 


PROGNOSTIC  DATA  REANALYSIS  ROUTINE 


UNITS 


SPEED  - MPS 

DIRECTION  - DEG 

HEIGHT  - METERS 

TEMPERATURE  - DEG  K (Virtual) 


RESULTS  FOR  LEVEL  L = ? (850  mb) 


LAT (OEG) 


t 


i 

f 


! 

DELETED  DATA 

I 

LON(DEG»  U<KTS)  V(KTS)  TEST  ! 


STATION  DATA 


NUMBER 

D VALUE 

temperature 

U COMP 

V COMP 

1 

29.00 

288.86 

10.99 

-.03 

2 

15.00 

283.76 

10.41 

.18 

3 

40.00 

290. 36 

8.2  3 

-2,87 

4 

-999.90 

-999.90 

-999.90 

-999.90 

5 

16.00 

289,26 

5.52 

.45 

6 

-999.90 

-999.90 

-999.90 

-999,90 

7 

-999. 90 

-999. 90 

-99T.90 

-999,90 

8 

-27.00 

289.96 

4.02 

-.01 

9 

29.00 

290.86 

1.24 

• 86 

10 

28.00 

288.  16 

4.39 

?«88 

~999.90  means  missing  data  from  Artillery  Met  Sections 
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GMC  DATA  ON  MINI  GRID* 


O VALUES 


AO 

AO 

39 

39 

39 

39 

39 

AO 

AO 

AO 

AO 

AO 

39 

39 

A1 

A1 

Al 

AO 

AO 

AO 

AO 

A2 

A1 

At 

Al 

Al 

Al 

Al 

A2 

A2 

A2 

A2 

A2 

At 

Al 

A3 

A3 

A2 

A2 

A2 

A2 

A2 

A3 

A3 

A3 

A3 

A3 

A3 

A2 

TEMPERATURE 


288 

289 

289 

289 

289 

289 

289 

289 

289 

289 

289 

289 

289 

289 

2 89 

289 

289 

289 

2 89 

289 

289 

289 

289 

289 

289 

289 

289 

289 

2 89 

289 

289 

289 

289 

2 89 

290 

289 

289 

289 

289 

289 

290 

290 

289 

2 89 

289 

289 

290 

290 

290 

U COMP 


1111112 
111112  2 

V COMP 

5 5 5 5 5 5 5 

5 5 5 5 5 5 5 

5 5 5 5 5 5 5 

5 5 5 5 5 5 5 

5 5 5 5 5 5 5 

5 5 5 5 5 5 5 

5 5 5 5 5 5 5 

* 

For  1230  GMT,  19  November  1974  as  interpolated  from  original 
^ GWC  prognoses  data. 

t 

i 

I 

I 
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* 

UPDATED  GWC  DATA  ON  MINI  GRID 


D 

VALUES 

P.O 

21 

36 

39 

41 

39 

32 

1 9 

24 

34 

32 

3 0 

27 

25 

29 

25 

29 

34 

33 

20 

22 

33 

29 

29 

30 

29 

25 

2 

36 

30 

29 

29 

26 

19 

7 

38 

32 

30 

27 

25 

20 

22 

35 

32 

37 

27 

26 

25 

24 

temperature 

289 

290 

291 

291 

290 

2 89 

288 

289 

290 

291 

291 

290 

2 89 

288 

290 

290 

291 

291 

290 

2 89 

283 

290 

290 

291 

290 

289 

2 89 

289 

290 

290 

291 

290 

268 

2 86 

285 

290 

290 

289 

289 

288 

286 

2 85 

289 

2 89 

289 

289 

288 

287 

287 

U COMP 

0 

-I 

-I 

-0 

1 

3 

4 

1 

1 

1 

1 

1 

1 

3 

2 

2 

2 

2 

2 

2 

3 

4 

4 

3 

4 

4 

4 

3 

6 

T 

5 

6 

7 

8 

8 

7 

7 

7 

6 

7 

8 

9 

6 

6 

3 

6 

7 

7 

6 

V COMP 

4 

3 

3 

1 

1 

1 

3 

4 

3 

2 

-0 

1 

0 

2 

3 

2 

2 

1 

2 

2 

2 

2 

2 

2 

1 

1 

1 

1 

0 

1 

3 

2 

1 

1 

1 

-1 

1 

1 

3 

2 

2 

2 

1 

1 

2 

2 

3 

3 

3 

* 

For  1230  GMT,  19  November  1974  after  modification  to  fit  Army 
observations  at  WSMR  taken  between  1210  and  1250  GMT. 


Appendix  H 


AAR  DATA  OUTPUT 


This  example  uses  data  from  updated  mlnlgrld  and  the  best  available  Army  rawlnsonde 
to  produce  a Computer  Met  Message  for  1230  GMT,  19  November  1974, and  for  the  location 
requested  in  the  input  parameters  (See  Appendix  F) . 


Ill 


f-  ^l»«  !■  ■ 

A«TILLeRY  APPLICATIONS  ROUTIN'? 


UNITS 

speeo  - MPs 
DIRECTION  - DECi 
HEIGHT  - METERS 

TEMPERATURE  - DEG  K (Virtual) 
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W 


'~v 


\ N 

300.0 
33S*0 

370.0 

405.0 

440.0 
4 75.0 
51  0.  0 

545 . 0 
5eo.  0 

615.0 
650.  0 

oB  j • 0 

7^0.  0 

755 . 0 
700.  0 
6P‘J.O 
060.  0 

505 . 0 

O 

065 . 0 

I 000. 0 


THt  HOH.nwlNL*  Gf?APH  ♦ = tJk-‘5T  A^AlLArtLr  SCUNJlNG  AM)  X = 5wC  lJO)AT-) 

X-INCPrMfNT  s P.OOO  AND  Y-lNfPr.MhNT  = -35,000 


- PO.  0 


•10.0 


0.0 


10.0 


, 0 


See  Figure  4 for  explanation  of  profiles. 
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★ 

5*  '.UN')  I N 


IN  THf.  r^LL'lAlNO  ♦ - U'^T  AVA|la‘L“  <'3U'J')T 

K-iNCPtUNT  = '->.000  AN'.)  Y- I NC^V  t. - N ^ 


lOvJ.O  ( 

J35.0  ( 

170.0  ( 

405.0  ( 

440*0  ( 

4 76.0  ( 

5 I 0 . O I 

546.0  ( 

580.0  ( 

015.0  ( 

650.0  ( 

085.0  ( 
7P0.0  < 

755.0  ( 

700.0  ( 

3PS.0  ( 

850.0  ( 

SO  1.0  ( 

O <0  *o  ( 

T • <)  ( 

1000.0  ( 


- 00.  > 


10*0 


o.  0 


IN  THt  rOLLTWING  GfvAPH  ♦ 


RFST  AVAILAjLf  SCUNOING  ANJ  X 


G«VC  u'^datlo  soundtn 


<- 

INCptm^-mT 

= 

000*000  AK^  Y- I 

PRFSS 

300*  0 

< 

335*0 

( 

370*0 

{ 

40S*0 

( 

440*0 

( 

4 75*0 

( 

♦ X 

5i0,0 

< 

4 X 

54  5 • 0 

( 

4 X 

seo.  0 

< 

•f 

X 

615*0 

( 

X 

650.0 

( 

♦ 

X 

605.0 

( 

720  • 0 

f 

♦ 

X 

755*  0 

( 

4-  X 

790.0 

< 

4 X 

825.0 

( 

«• 

X 

860.  0 

( 

X 

8 95.0 

( 

930*0 

( 

965*0 

( 

1 000*0 

( 

- .15  • 00  0 


I **O00.  0 
H(  I 


0 


0 


4000  *0 


0000*0 


I 2000*0 


STANOAHO  FORMAT 


PRESS 

1 N MHS 

« TEMP 

IN  DEGIK ) . 

U AND  V 

COMP 

IN  MPS 

BEST  AVAILABLI- 

SOUNDING* 

UPDATED  SOUNDING  ** 

- ZONE 

PRESS 

TEMP 

UCQMP 

VCOMP 

ZONE 

PRESS 

TEMP 

UC  OMP 

VC  OMP 

1 

876*0 

287.0 

3.6 

2.0 

1 

876.0 

289.0 

-3.7 

• J 

2 

866*  0 

286.6 

9.  8 

.3 

? 

H6<i.  0 

288.8 

3. 1 

-1  .4 

3 

840.0 

285.2 

12.3 

.4 

H40.  0 

288.  1 

H.O 

-1.1 

4 

801  • 0 

282.0 

15.3 

1 .8 

4 

801.0 

284.2 

12.4 

• 1 

5 

753.0 

278.7 

18.4 

-2.2 

« 

751.  0 

280.0 

1 7.3 

-1.9 

6 

708.0 

275.2 

12.8 

-5.3 

6 

70H.  0 

275.7 

13.4 

-7.0 

7 

666.0 

273.  1 

1 5.  4 

-1.5 

7 

666 . 0 

273.4 

15.3 

-3.2 

8 

625.0 

270.8 

1 8.  1 

4.0 

8 

625.0 

270.9 

16.9 

2.3 

F- 

587.  0 

268.9 

17.0 

.3 

9 

587.0 

268.9 

14.7 

- 1 .4 

i 10 

550.0 

267.  1 

21.5 

-2.5 

10 

550.0 

267.0 

18.2 

-4.2 

1 11 

516.0 

264.  1 

19.5 

-1.2 

1 1 

516.0 

263.  8 

14.9 

-2.9 

12 

484.0 

259.8 

1 8.4 

-4  .« 

12 

464.  0 

259.5 

13.8 

-6.5 

LATITUDE  LONGITUDE 

-105.23  (Midpoint  of  Trajectory) 

COMPUTER  MET  MESSAGE  FORMAT 

PRESS  IN  MRS.  TEMP  IN  TENTHS  OF  DEGIKI.  DIRECTION  IN  TENS  OF  MILS.  S"EED  IN  KNOTS 


BEST  AVAILABLE  SOUNDING  UPDATED  SOUNDING 


ZONE 

PRESS 

TEMP 

OIR 

SPEED 

ZONE 

PRESS 

TEMP 

OIR 

SPEED 

0 

8 76 

2870 

427 

8 

0 

876 

2890 

169 

7 

1 

866 

2866 

477 

19 

1 

866 

2868 

5 24 

7 

2 

64  0 

2852 

477 

24 

2 

840 

2881 

497 

16 

3 

801 

2820 

468 

30 

3 

801 

2842 

4 79 

24 

4 

753 

2787 

492 

36 

4 

753 

2800 

50  2 

34 

5 

708 

2752 

520 

27 

5 

708 

2757 

529 

29 

6 

666 

2731 

490 

30 

6 

666 

2734 

501 

30 

7 

625 

2708 

458 

36 

7 

625 

2709 

466 

33 

8 

587 

2680 

478 

33 

e 

587 

2689 

490 

29 

9 

550 

2671 

492 

42 

9 

550 

2670 

50  3 

36 

10 

516 

264  1 

486 

38 

10 

51  6 

2639 

4 39 

29 

11 

4 64 

2598 

506 

37 

1 1 

484 

2595 

525 

30 

northing®  3. 1 2T48A ASOFOe  EASTING®  1 . 5 45 3 3 7 770* 07 
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